1 



Confluent and Double-Confluent Heun Equations: Convergence of Solutions in 

Series of Coulomb Wavefunctions 

Lea Jaccoud El-JaiclJl] and Bartolomeu D. B. FigueiredcH 
Centro Brasileiro de Pesquisas Ffsicas (CBPF) 
Rua Dr. Xavier Sigaud, 150 - 22290-180 - Rio de Janeiro, RJ, Brasil 



(N 

o 

<D 
CO 

o 

(N 



Abstract 

The Leaver solutions in series of Coulomb wave functions for the confluent Heun equa- 
tion (CHE) are given by two-sided infinite series, that is, by series where the sum- 
mation index n runs from minus to plus infinity [E. W. Leaver, J. Math. Phys. 27, 
1238 (1986)]. First we show that, in contrast to the D'Alembert test, under certain 
conditions the Raabe test assures that the domains of convergence of these solutions 
include an additional singular point. Further, by using a limit proposed by Leaver, we 
obtain solutions for the double-confluent Heun equation (DCHE). In addition, we get 
solutions for the so-called Whittaker-Ince limit of the CHE and DCHE. For these four 
equations, new solutions are generated by transformations of variables. In the second 
place, for each of the above equations we obtain one-sided series solutions (n > 0) by 
truncating on the left the two-sided series. Finally we discuss the time dependence 
of the Klein-Gordon equation in two cosmological models and the spatial dependence 
of the Schrodinger equation to a family of quasi-exactly solvable potentials. For a 
subfamily of these potentials we obtain infinite-series solutions which converge and are 
bounded for all values of the independent variable, in opposition to a common belief. 
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1 Introduction 

In 1986 Leaver |T] found two types of solutions in series of confluent hypergeometric 
functions for the confluent Heun equation (CHE) and presented a limit procedure to generate 
solutions for the double- confluent Heun equation (DCHE) out of solutions for the CHE. 
Later on [21 EJ H] we have found that there are two other physically relevant equations whose 
solutions can also be derived from the Leaver solutions for the CHE and DCHE by means of 
a procedure called Whittaker-Ince limit. Further, from solutions of the CHE and/or DCHE, 
we can find solutions for the Mathieu, Whittaker-Hill and spheroidal equations [lj Ej. 

In view of the above connections, from the inception we establish the convergence prop- 
erties of Leaver's solutions. We consider only the expansions in series of Coulomb wave 
functions which are given by a set of three solutions, one in series of regular confluent hy- 
pergeometric functions and two in series of irregular functions. By redefining the Coulomb 
functions, we avoid difficulties arising from the Leaver definitions and find that the con- 
vergence of the solutions for the CHE and its Whittaker-Ince limit follows from the Raabe 
test. We investigate as well the transformations of the Leaver solutions and their limiting 
and special cases. 

First we write the aforementioned equations, present the connections among them and 
call attention for the fact that there are three types of series expansions whose convergence 
require different treatments. After that we introduce the D'Alembert and Raabe tests for 
convergence and outline the structure of the article. 

The Leaver form for the CHE is [1] 

d 2 U dU 
z(z - zo)-^ + (#i + B ^ z )~^ + l B 3 ~ ^{z - zq) + co 2 z{z - z )\ U = 0, (1) 

where Bi, r\ and u are constants. The points z = and z = zq (if z ^ 0) are regular 
singular points, whereas z = 00 is an irregular point. Since z is free, by taking z = 
Leaver obtained the DCHE 



z 



d 2 U , „ „ x dU 



(B 1 + B 2 z) — +(B 3 - 2 V uz + u 2 z 2 ) U = 0, [B x ^ 0, u ^ 0] (2) 



dz 2 dz 
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in which case z = and z = oo are both irregular singularities. 

In addition, the CHE and the DCHE admit a (Whittaker-Ince) limit which changes the 
nature of the singularity at z = oo, keeping unaltered the other singular points. This limit 
is given by [21 [3] 



uj — > 0, 77—7- 00, such that 2t]u = —q, (Whittaker-Ince limit) 
where q is a nonvanishing constant. The Whittaker-Ince limit of the CHE is 

z{z -*o)^ + (B 1 + B 2 z)^- + [B 3 + q(z - z )} U = 0, (g ^ 0) 



(3) 



(4) 



(if q = this equation can be transformed into a hypergeometric equation), while the limit 
of the DCHE is 



z 71 + ( 5 i + 5 2^)— + S3 + g* c/ 

dz z dz 



0, (g ^ 0, B 1 ^ 0) 



(5) 



(if g = and/or B\ — the equation degenerates into a confluent hypergeometric equation 
or simpler equations). This last equation also follows from Eq. (J4j) when z = (Leaver's 
limit). 

The Mathieu, Whittaker-Hill and spheroidal equations have been studied by themselves, 
but they are particular cases of the above equations. The Mathieu equation reads [B] 



d 2 w 
du 2 



a 2 [a — 2k 2 cos(2au)]w = 0, (Mathieu equation) 



(6) 



where a and k are constants, while a = 1 or a = i for the Mathieu or modified Mathieu 
equation, respectively. This equation is transformed into particular instances of Eqs. ([T|), 
(j2D and (J4j) by the substitutions of variables given in Eqs. ( 1631) . f llOOl) and ( 180|) . respectively. 



(7) 



The Whittaker-Hill equation (WHE) can be written in the form [7J |8] 



d 2 W 



du 2 



r) 



< 2 ~ (p + 1)£ cos(2^u) + cos(4?tt) 



W = 0, (WHE). 



where i9, ^ and p are parameters; if u is a real variable, this represents the WHE when c = 1 
and the modified WHE when q = i. The WHE reduces to the CHE (fTj) and DCHE fl2J by 
the substitutions (1511) and (j9oT) . respectively. Finally, the spheroidal equation reads [2j 



d_ 

dy 



(1 - v 2 ) 



dS{y) 
dy 



+ 



A + 7 2 (l-y 2 ) 



5(2/) = 0, 



where 7, A and /i are constants. The substitutions (158aj) transform this into a special case 
of the CHE (UJ. 

On the other side, in general there are solutions given by three different types of se- 
ries, called two-sided infinite series, one-sided infinite series and finite series. These take, 
respectively, the forms 



N 



9n{Z) 



(9) 



n=0 



n=0 



where f^{z) and g n {z) are functions of the independent variable z, N is a non-negative 
integer and v is a parameter which does not appear in the differential equations. In the 
present case, the series coefficients a n and b n satisfy three-term recurrence relations. No 
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finite-series solutions are known for the Mathieu equation (|6]) nor for the Whittaker-Ince 
limit (ED of the DCHE. 

Two-sided infinite series, the only considered by Leaver, are necessary to assure the 
convergence of solutions of equations in which there is no free parameter, as in scattering 
problems [2] or in some wave equations in curved spacetimes [10] . In such cases the parameter 
v must be determined as solutions of a transcendental (characteristic) equation. However, 
when truncated on the left (n > 0), the two-sided infinite series give one-sided infinite series 
which are useful for equations having a free parameter; in turn, these lead to solutions given 
by finite series if the parameters of the equation satisfy certain constraints. 

Finite-series solutions are suitable for quasi-exactly solvable (QES) problems, that is, 
for quantum-mechanical problems where one part of energy spectrum and the respective 
eigenfunctions can be computed explicitly [TT| IT2] . For QES problems obeying equations of 
the Heun family [3], that part of the spectrum may be derived from finite-series solutions if 
these are known. Indeed, a problem is also said to be QES if it admits solutions given by finite 
series whose coefficients necessarily satisfy three-term or higher order recurrence relations 
[T3] . and is said to be exactly solvable if its solutions can be expressed by hypergeometric 
functions. 

The convergence of two-sided infinity series is obtained from the limits 



By the D'Alembert test the series converges in the intersection of the regions for which 
L\ < 1 and L 2 < 1, and diverges otherwise (if Li = L 2 = 1, the test is inconclusive). 
Leaver's definitions for the Coulomb wave functions lead to ratios between terms presenting 
square roots (except if rj = 0) which make difficult to deal with the convergence tests. To 
avoid this problem, we use alternative definitions that, in addition, permit to apply the 
Raabe test for the solutions of the CHEs. By the Raabe test [T4"l [T5j . if 



where A and B are constants, then the series converges in the region where A < — 1 and 
B < — 1, and diverges otherwise (if A = B = — 1, the test is inconclusive). For one-sided 
series the convergence may be enhanced since we use only the limit L±, while for finite series 
the convergence must be decided from the behavior of each term of the series. 

Furthermore, by using transformations of variables we find four sets of two-sided solutions 
instead of one set as in Leaver. By the Raabe test, under certain conditions these solutions 
converge absolutely for \z\ > \zq\ or \z — z Q \ > |z | rather than for \z\ > \zq\ or \z — z \ > \zq\; 
the one-sided solutions given by series of regular confluent hypergeometric functions converge 
for \z\ > 0. Nevertheless, the behavior of each solution for z — > oo must be analyzed carefully 
because, in computing L\(z) and L 2 (z), we assume that z is bounded. We have also to 
examine the behavior of the solutions at the finite singular points because the series appear 
multiplied by factors which may become unbounded at such points. 

Taking into account the above considerations, we proceed as follows. In Section 2 we deal 
with the two-sided infinite expansions for CHE ([I]), its limiting and particular cases. For the 
spheroidal equation, we get the Meixner solutions in series of Bessel functions [Sj instead of 
the Chu and Stratton solutions [16] mentioned by Leaver. 

In Section 3 we apply the analysis of section 2 to one-sided expansions. These are 
generated by truncating on the left-hand side the two-sided series: by requiring that n > 
we determine the parameter v of the two-sided solutions as function of the parameters of 
the differential equations. For r/ ^ we find three types of recurrence relations for the series 
coefficients; for 77 = 0, only two types. Conditions for finite-series solutions are obtained. 




(10) 
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In section 4, we consider examples which illustrate the consequences of Raabe test. In 
particular, we show that the Schrodinger equation for some quasi-exactly solvable poten- 
tials admits infinite-series solutions which are convergent and bounded for all values of the 
independent variable. Thus, in addition to the energy levels resulting from finite series, in 
principle it is possible to get additional energy levels as solutions of characteristic equations 
corresponding to the infinite series. 

In section 5 there are some final considerations and remarks on open problems concerning 
the Leaver solutions. Appendix A gives the normalization used for the Coulomb functions 
and takes the case t] = as a criterion to decide in favor of one of two possibilities for the 
ratio between successive Coulomb functions. The derivation of the recurrence relations for 
one-sided and two-sided series is given in Appendix B. Appendix C lists one-sided solutions 
for the CHE, omitted in section III. 



2 Two-sided Series Expansions 

In this section we examine separately the two-sided expansions for each of the four equa- 
tions: CHE, DCHE and their Whittaker-Ince limits. In the Whittaker-Ince limit, the ex- 
pansions in series of Coulomb functions give solutions in series of Bessel functions, excepting 
the solutions fllQ4al) for the Whittaker-Ince limit of the DCHE. The solutions of the CHE 
with rj = are also expressible by series of Bessel functions. In all cases, given an initial 
set of solutions, new sets are generated by transformations of variables which preserve the 
form of the differential equations. Notice that the linear independence of the functions used 
as basis for the series expansions will impose restrictions on the characteristic parameter v 
and/or on some parameters of the differential equations. 

In Eqs. ( 16011621) we recover the Meixner solutions for the spheroidal equation as particular 
cases of the solutions for the CHE, while in Eqs. ()81aH81cj) we recover the usual solutions 
in series of Bessel functions for the Mathieu equation as particular cases of the solutions for 
the Whittaker-Ince limit of the CHE. In both cases the Raabe test modifies the domains of 
convergence of some solutions. 



2.1 Solutions for the CHE 

The initial set of solutions, Ui(z), is reconstructed in Appendix B. It reads 

UiO) = z~^ b l^n+u (v, wz) , W n+V (rj, uz) = (<f) n+u , ip+ +v , ip~ +u ) (rj, uz) , (12) 

n 



where denotes two-sided series, 4> n + v and ip n+v represent the Coulomb wave functions 
defined in Eqs. (1A.9I) and ( 1A.10I) . and the coefficients satisfy three-term recurrence rela- 
tions. So, we have a set of three expansions, one in series of regular confluent hypergeometric 
functions and two in series of irregular functions. This set corresponds to Leaver's solutions 
who have used the definitions ( 1A.11I) and ( 1A.12[) for the Coulomb functions. In addition, if 



U (z) = U(B±, B 2 , B 3 ; z , cu, r/; z) denotes an arbitrary solution of the CHE, we can find other 
solutions by means of the transformations T\, T 2 , T 3 and T4 which operate as [3j |5] 

T x U(z) = z 1+B ^ z °U(C u C 2 , C 3 ; z , u, rj; z), 

T 2 U(z) = (z-z ) 1 - B *- B i/ z °U(B 1 ,D 2 ,D 3 ;z ,u,r)-,z), 

v J 

T 3 U(z) = U(B 1 , B 2 , B 3 ; z , -u, -77; z), 

T 4 U(z) = U(-B 1 - B 2 z , B 2 , B 3 + 2t)uZq, z , -u,m z - z), 
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where 

C l = -B 1 -2z , C 2 = 2 + B 2 + ^ C 3 = 5 3 +(l + f)(fi 2 + f) 



D 2 = 2 - B 2 - D 3 = B 3 + ^(^ + B 2 -l 



(14) 



These transformations allow constructing a group with 4 sets of two-sided series U, (i = 
1, • • ■ ,4) where coefficients 6^ satisfy recurrence relations having the form 

<K +1 + 0% + = 0, [-00 < n< oo] (15) 

where a l n , /?* and 7^ depend on the parameters of the differential equation as well as on v and 
n. These relations lead to transcendental (characteristic) equations given as a sum of two 
infinite continued fractions. By omitting the superscripts of a l n , f3 l n and 7^, the characteristic 
equations read 

_ a_i7o a_27-i "-37-2 «o 7 i "i 7 2 «273 , . 

^^I^IP" 4 (16) 
which are equivalent to the vanishing of the determinants of infinite tridiagonal matrices, as 
in Eq. (1321) . If the CHE has no free parameter, Eq. (Tl6l) may be used to find the possible 
values of v (characteristic parameter); if the CHE has an arbitrary parameter, Eq. ( fl6|) 
permits to find the values of that parameter corresponding to suitable values of v. 

To analyze the properties of the solutions we write explicitly each of the three solutions, 
instead of using the abbreviated form (1T2"|) . Thus, we denote by Ui = (Ui, , C/jf) the 
solutions associated respectively with the functions (4> n+u , i^n+vi VWi/)- This gives the 
solutions f fT9|) which, by the transformations f JT3|) . generate the three sets of solutions that 
have not been considered by Leaver. The four sets of two-sided solutions are denoted by 

Ui(*)= [Ui(z),Uf(z),Ur(z)] , z = l,-..,4, (17) 

if rj 7^ 0; if 77 = the notation is given in Eq. (J47J). From these solutions we get eight sets of 
one-sided solutions which are denoted by 



Ui(z),U?(z),Ur( z ) , i = l,... ,8, (18) 



and do not depend on v. In fact they are generated by expressing the parameter v of each 
as two different functions of the parameters of the CHE. The convergence of the solutions U; 
is obtained by considering only the limits n — > 00 in the computations given in this section. 

2.1.1 The four sets of the solutions 

Explicitly the first set Ui, given in Eq. f fT2|) . reads 



b 1 \2iui z\ n+u+1 ~ 

U x (z) = e luz ^pr " t-$ [n + u + l+ ir],2n + 2u + 2; -2iuz\ 

.T I 2 v I —\ 



U±( z ) = e ±iuiz V i — \n + v + 1 ± iri, 2n + 2i/ + 2; ^2iu;zl 

n L J 

where, in the recurrence relations ffT5]) for 

2iaj2o \n+v+2— ^2. 

a 1 



(19) 



(2n+2v+2) (2n+2i/+3) 



ft = -a, - W2o - („ + „ + 1 - f ) („ + „ + f ) - 9g=gfe|l , (20) 

1 2iwzo [n+i/+^- — l] + (n+v+ir?)(n+i/-«?7) 

%i ~~ (2n+2i/-l)(2n+2i^) ' 
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By applying the transformation T3 on Ui, we find the equivalence 

T 3 [U 1 (z),U+(z),U{(z)] [U l {z)^{z),Ut{z)]. (21) 

Precisely, we find T 3 [U 1: U? t C/f] = [t/i, f/f, t^] with 

U x {z) = e™ z n r r 2n + 2 z/ + 2l $ (™ + ^ + 1 + ^ 2 ™ + 2 ^ + 2 ! ~ 2 ^ z ) > 

E^l f2zo;zl n+1/+1 ~ B2 / 2 
=-tf (n + v + l±iri,2n + 2v + 2; ^-2iuz) , 
r \n + z/ + 1 =F tri] 

where the recurrence relations for b\ are 

Up to a multiplicative constant independent of n, we can set 6* = {—l) n h l n in order to 
establish relation ( 12TT) . Thus, the transformation T3 is ineffective in the present case. The 
remaining transformations allow to form a group constituted by four sets of solutions, namely, 

V^z), U 2 (z) = T 2 V 1 (z); U 3 (z) = TJJ^z), U 4 (z) = T 4 U 2 (z) = TJJfa). (22) 

In fact, from the explicit forms of these sets, one verifies that T\ does not generate new 
solutions when applied on Ui and U2; similarly, T 2 has no effect on U3 and U4. 

The three sets of solutions generated by the preceding transformations are written below. 

b 2 \iujz\ n 

U 2 (z) = f 2 (z)e^ z Yl r[ +2z/ + 2] $ [n + V + 1 + iv ' 2U + 2U + 2; ~ 21UZ] 

n (23) 
h 2 \—2ii )7~\ n 

where / 2 = f 2 (z) = ^+(BiAo)+(s 2 /2) ^ _ z ^i-b 2 -(Bi/ Zo )_ The coefficients for the recurrence 
relations are given by 



2^ [n. + , + ^][n + , + l + ^ + ^] 

(2n+2i/+2)(2n+2i/+3) ' ^« 

Br a ( 24 ) 

~~ (2n+2i/-l)(2ri+2f) " 

The transformation T 4 acting on Ui gives the set U3 = 113(2;), namely, 

U 3 = / 8 e*"[*-* »] £ r[2r+^+2l W $ [n + " + 1 + ^ 2n + 2 " + 2; ~ 2 ^ " Zo)] ' 

n 

vf = he ±M>-M y ft CM*-*)]- (25) 

* [n + 1/ + 1 ± irj, 2n + 2i/ + 2; =F2iw(z - z )] , 
where f% = fs{z) — (z — zo) u+1 ~^ B2 ^ 2 ^ and, in the recurrence relations, 



2^ [n + , + 2-^][n+, + l+^ + ^] 
"n (2n+2f+2)(2n+2f+3) ' i J m 

g 2iu)zo[n+u— 1+-^-] [n+y— ^ — -j^-J (n+i/+i77)(n+J>— £77) 
Tra ~~ (2n+2i/-l)(2n+2i/) ■ 



(26) 
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The fourth set, obtained by applying T\ on U3, reads 

U A = h e^-M rfo^lV^oi" * [" + " + 1 + ir >> 2n + 2 » + 2 ; ~ *>)] , 

1 2n + 2z/ + 2 

n J 

^ =/4e ^-, IE jfPMf^£)r x (27) 

\l> [to + z/ + 1 ± irj, 2n + 2v + 2; =f2zw(z - z )] , 
where ^4 = fi(z) = z l+<yBl ^ Z0 ^ (z — zq)""^ 2 ^ 2 ^^ 1 ^^ and, in recurrence relations for 6^, 



1 



2iwzo [ra+iH-^p-] 



6 A = B l 



n (2n+2u+2){2n+2u+3) ' y n l^ni 

4 2itJ2o [n+l>+l— ^] [n+H-^ + ^J (n+iH-i»7)(n+i/— 117) 
~~ (2n+2u-l){2n+2v) ' 

If there is no free parameter in the CHE, v must be determined as solutions of a charac- 
teristic equation. However, by considering the form of the solutions and respective recurrence 
relations for the series coefficients, we find that 

2v and v =F it] cannot be integers (29) 

for two-sided series. The restriction v^ir] ^integer assures that factors 1/T(to + v + 1 ± irj) 
which appear in U^(z) are not zero for any value of n; assures as well that the factors 
(to + v + irj) (n + v — irj) in 7^ do not vanish for any n. In fact, such restriction is necessary 
to have two-sided infinite series for the three solutions in each of the four sets Uj. 

The condition 2v ^integer is necessary in order to avoid two terms linearly dependent 
in the series of Uf{z). Indeed, suppose that 2u =integer in the solutions U^(z). These are 
series expansions in terms of 

B±(z) = [-2iojz] n+v+l ^[n + u + l±ir],2n + 2u + 2; ^2iuz]. 

By setting to — ri\ and using (1A.3I) . we find 

£*(z) = ±(-\y TU [-2iuJz\~ ni ~ u ^[-n 1 -u±ir], -2n x - 2u; =p2iwz]. 

Hence, B„ and B^ 2 are proportional to each other for some to = n 2 such that TOx + n 2 + 1 = 
— 2v. Similar results are found for the other solutions \jf[z). On the other side, by supposing 
that 2v ^integer, the functions $(a, c; y) which appear in Ui(z) are well defined because the 
parameter c = 2n + 2v + 2 cannot be a negative integer. Nevertheless, see in the paragraph 
containing Eqs. f )55|) and f )56|) some remarks concerning the case rj = 0. 

According to Eqs. ( IA.4j) . if the conditions ( 129]) are true, the three hypergeometric func- 
tions are linearly independent and each one can be written as a combination of the others 
by means of ()A.5j) . In this case, in a common region of validity, we can write one solution 
of a given set as a superposition of the others. However, the three series are really doubly 
infinite (—00 < to < 00) if, in addition to f )29|) . v satisfies the restrictions 

v ± ^ and v ± + are not integers. (30) 

These conditions assure that a l n and Y n do not vanish for any value of to. In effect, if a % n = 
for some to = N\, the series should begin at to = N% + 1 in order to assure the validity of the 
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theory of the three-term recurrence relations; by the same reason, if 7^ = for some n = N2, 
the series should terminate at n = — 1. 
Notice that, for two-side solutions, 



1 

n ' 



a n r Yn+l 



a nln+li 



2,3,4]. 



(31) 



Thence, Eq. (fl6|) implies that all the solutions satisfy the same characteristic equation 
and, consequently, the parameter v takes the same values in all solutions. In addition, as 
noticed by Leaver, the characteristic equations are periodic in v with period 1. In effect, 
in order to indicate that the coefficients depend on v we rewrite the recurrence relations as 
&nK+i + PnKi + IrMi-i = or as the following tridiagonal matrix equation: 



In In "n 

ryV RV n V 

ln+1 f^n+l "n+l 

7n+2 Pn+2 a n+2 



If 
ii 

b" 

n 



n+1 







[—00 < n < oo] 



(32) 



where denotes the null column vector. Thence, the values for v may be determined by 
requiring that the determinant of the above matrix vanishes. However, as 



In i n+1' rn r'n 



u+1 



+1' 



<y u+1 = a v , 



and — oo < n < oo, the matrix and its determinant are not modified by the replacement 
v — > v + 1 (or v — > v + N, where N is any integer) . 

Further, if (|29~I) and (130]) are fulfilled, all coefficients can be written in terms of b\. Up to 
multiplicative constants independent of n, we have 



r[n+u+^} r[n+,+i+a 



1 _i_i!2 

2 



b 3 



(_!)» r[n+,+l-^-^] ^ 



4 _ (-1)" r[n+,+2-^] x 



(33) 



As an example, we consider the solutions W(u) for the WHE ([7]). These may be obtained 
from the solutions U (z) of the CHE (pQ) by taking 



W(u) = U(z), z = cos 2 (qu), [q = 1, i] 
R - 1 ft — 1 R - (g+jOizjj , 



2' 



17] 



p+1 
2 



WHE as 
a CHE. 



(34) 



Thus, the solutions = (Ui, U^, U i ) lead to four sets of solutions 

Wi(ii) = [Wi(«) = E^*), = ^(z), = Ur{z)} . (35a) 

In this case, the coefficients of the recurrence relations for b\ simplify to 

«n = f, n = -B 3 -r)cu- [n + u+\]\ 1 l = -ifln + u + ir 1 \[n + u-ir,l (35b) 
whereas Eqs. (I33|) reduce to 



b 2 n = (n + v + \) b\, 



b 3 
n 



(_l)n (n + ^+|) 



(35c) 
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2.1.2 Convergence and an asymptotic behaviour 

The D'Alembert test implies two subgroups of solutions since Ui and U2 converge for 
any finite z such that \z\ > \z \, whereas U3 and U4 converge for \z — Zq\ > \zq\. However, by 
the Raabe test they may converge also at \z\ = \zq\ and \z — zq\ = \zq\ under the conditions 





Re 


B 2 


+ 


§1 
z ° . 


\z\ > \zq\ if < 










Re 


B 2 


+ 


§1 
z _ 



< 1 in Ui, 
> 1 in U 2 ; 



Zo\ > \ z o\ if 




^ > -1 in U 



3 j 



< -1 in U 



(36) 



■1; 



where the restrictions on parameters of the equation are necessary only to assure convergence 
at \z\ = \zq\ or \z — zq\ = \zq\. In particular, for the solutions of the WHEs we find 



cos(cu)| > 1 in Wi, 
sinfcw)! > 1 in W3, 



I cos(?m)| > 1 in W 2 , 
sin(«j«)| > 1 in W4. 



(37) 



Thence the two-sided solutions are useless for the WHE (that is, for q = 1 and w=real), 
but may be useful for the modified WHE (q = i, w=real). If Re[_B 2 + (Bi/z )] = 1 and 
Re[-Bi/z ] = — 1 in (1361) . the Raabe test becomes inconclusive in the sense that the solutions 
may converge or diverge at \z\ = \zq\ or | z — zq \ = \zq\. 

To obtain the conditions (1361) it is sufficient to consider the convergence of the first set of 
solutions. The results for the other sets arise from transformations ffTB"]) applied in the order 
given in 022)) . Thus, by using the form (]12p for the first set, the domains of convergence 
follow from the ratios 



lim 

n— >oo 



b l n+1 W n+v+1 {r],uz) 



b x n 'W n+v {'q,uz) 



and 



lim 

n— >— 00 



6i_ 1 ^ n+I/ _ 1 (77,W«) 



The ratios b\ +1 /b\ and b\ l _ x jb\ come from the relations + + 7„&„_ 1 = which, 

when n — > ±00, yield 



UjJZq 
—IUIZqU 



+ f -|)+0(^)]%i-2n [n + 2u + 1 + 0(1)] 



n 



2v + B 2 + ^-\ + 0(±) 

z zn 2 V n I 



bl 



Hence, the minimal solution for b\, x jb\ when n — > 00 is 



"n+l UZQ 

bl 2i 



and the minimal solution for bl l _ 1 /bl l when n — > —00 is 



0. 



(3* 



(39a) 



bl 2n 2 



b j+i 2n 

! — 



bl 



■2 

VjJZQ 



1 + M 2v + B 2 



Si 

zo 



(39b) 



On the other hand, from relations (1A. 16|) and (|A.17[) we find that, for finite values of . 



n — y 00 : 
n — > —00 : 



<Pn+v 



— 71 



~ [1 " f 1 

o;z L in J 



(40) 



Thence, by means of (I39ap . we find 
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and, by means of ( I39bl) . 



n — > — oo 



Bo 



--o 



(0 n+ „Vn+,)- (41b) 



From these limits we get 



lim 

n— ¥00 



n +l</>n+i/+l 



U1 2 ZqZ 



An 2 



lim 



Pol 
\z\ 



1 + -^Re ( B 2 - 2 + — 



and 



lim 

n— >oo 



± 



lim 



\Zq\ 
\z\ 



1 + ^Re ( B 2 - 2 + ^ 
\n\ V 



(42) 



(43) 



So, by the D'Alembert test the series converge absolutely for \z\ > \zq\ because the right- 
hand sides of ( 142]) and (143"]) are < 1. However, if |z| = \z \, by the expressions ( ITT]) for the 
Raabe test, the series converge even for \z\ = \zq\ provided that the numerators of \n\ in ( 142]) 
and (1431 are < — 1, that is, 

if Re [B 2 + (Bi/zo)] < 1, the series in Vi(z) converge for \z\ > \z \. 

If Re [B 2 + (Bi/z Q )} > 1, the series diverge and, if Re [B 2 + (Bi/z )] = 1, the test is in- 
conclusive. The convergence regions f[3"6"]) for the other sets of solutions are obtained by 
transforming the parameters and the variable z of Ui in accordance with Eqs. (1221 . Only 
the limit n — > oo is relevant for one-sided series (n > 0) and, then, the solutions Ui converge 
for any finite value z in virtue of the first limit given in (142]) . The convergence of Iff is 
similar to that of Uf. 

Since the previous regions of convergence were derived by supposing that z is finite, now 
we consider the behaviour of the solutions at z = oo. By using (1A.7D we find that, when 
z — > oo, 



Ui(z) 



-IT]- 



E 



r[n+i/+l~ irj\ 



r i ■ Bl * — * ( -\\n,— v—\-\-vniX 

L J ^— ' T[n+v+l+ir)\ 



(44) 



Thus, U\ may be unbounded by virtue of the exponential factors. This is consistent with 
the fact that, if conditions (129]) and ( 130]) are satisfied, then U\(z) can be written as a linear 
combination of U±(z) and U^(z). In fact, when z — > oo, Eq. (IA.6j) gives 



Uf(z) 



-2iuz] 



-17]- 



T[n+v+l—ir)\ ' 



< arg(-2^) < f , 



So / i \n~-is — l-f-ZTi J-.1 

ITT W ~ [-2iuz]<^ + £ ( "r[U^i-H^ 



^ < a,ig(2iujz) < ^. 



(45) 



Thus, the series in converge at z = oo but one of them may be unbounded depending 
on the exponential factors. For instance, if Re(iuz) — > oo, U± — > oo but f/f is bounded. 



2.1.3 The case rj — 0, the spheroidal and Mathieu equations 

Taking 77 = and keeping fixed the other parameters, the previous solutions are rewritten 
in series of Bessel functions of the first kind, J K (y), and in series of the first and the second 
Hankel functions, H k(u) and H^(y). We also include the Bessel functions Y R (y) of the 
second kind. These four functions are denoted by Z n{y) - or by ^^(y) - according as 

ZW(y) = J K (y), Z&>(y)=Y K (y), Z^(y)=H^(y), Z^(y) = H^(y). (46) 
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There are connections among these functions. For example, the relation Y K = [H^ — 

(2,) 

H v]/(2i) permits to obtain the expansion in series of Y K as a linear combination of the 
expansions in series of Hankel functions. Thus, we get four sets of solutions, each containing 
four solutions. These sets are written as 



UH\*) = [U (1 Az),U^(z),Ur(z),U ( :>(z)\ , i = 1,2,3,4, (47) 

where the right-hand side corresponds to the Bessel functions f|46|) . For one-sided series there 
are eight sets U^. The solutions Ui lead to U± which, in turn, give the other sets by means 
of the transformations (1221). that is, 



r (l)/ X rr(2)/ x rr(3), x j Ai) 



U%\z) = T 2 U^(z); U%\z) = TtU^iz), U^{z) = IW^z). (4* 
Thus, we put ?7 = in Ui given in ( fl9l) . use the relations ( 1A.19j) together with [17] 



r(2z) = 2 2 ^ 1 r(z)r[z + (i/2)]/ v ^, 

and redefine the coefficients as a* = 2 n 6*/r(n + z/ + 1). So, we find 
In the recurrence relations f|T5|) for a*, we have 

! ^,[»+^24][n+,+l-|4] 
fl n (2n+2i/+3) ' 

/^ = -(n + z/ + l-f) (n + r/+f) -5 3 , (49b) 

1 **o[n+*+^-l] ^ + 

^ri (2n+2i/-l) 

The other sets are given by (/3* = 2z> 7^ 0, ±1, ±2, ■ ■ • ): 



= Z^ + ^[z - Zot B ^J2 a n Z nL^ 

n 

2 _ ^o["+^]["+^+ 1 +^ + ^] 2 _ ^o[n+^+l-^][n+^-^- 

a n (2n+2v+3) ' (2»+2i>-l) 



a n (2n+2i/+3) ' %i (2n+2i>-l) 



3 w^ol 



4 _ ^»K^+^]h+^+l-^-^] 4 _ ^ [n + ,+l-^][n + ,+ |L + ^] 

a n (2n+2^+3) ' In (2n+2u-l) 

By using the relations [T7] 

J K {ye™) = e™ K J K (y), Y K {ye™) = e~^ nK Y K (y) + 2i cos(™) J K {y), 

H { 1 ] {ye m ) = -e-^H^iy), H (2 J (ye™) = e™ K H {1 ^(y) +2cos(irK)H®(y) 



(50) 



(51) 



(52) 



(53) 
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with k = n + v + (1/2), we find that the change u — > —u does not lead to new independent 
solutions. In this sense, once more the transformation T 3 is ineffective. 
Also in the present case (77 = 0) conditions (130|) . that is, 

v ± ^ and v ± (j^ + are not integers, [ see relations (130]) ] (54a) 

are necessary in order to have two-sided infinite series, and relations ( 1331) hold for the 
coefficients a\. On the other side, the restrictions (1291) are replaced by conditions 2v 7^ 
0, ±1, ±2 • • • which assure the independence of the Bessel function: 

2u 7^ 0, ±1, ±2, • • • , [independence of Bessel functions]. (54b) 

In fact, it is necessary that v 7^ ±1/2, ±3/2, • • • in order to avoid two linearly dependent 
functions of integer order, like Zg\y) and Z^J(y) [Zg = (— VfZ_H\, where t is zero or positive 
integer. In addition, v 7^ 0, ±1, ±2, • • ■ assures the independence of the Hankel functions in 
the same series: on the contrary, we would have functions like H^ + /ii 2 ) and H_£_(i/ 2 ) which 
are proportional to each other since [IB] 



tf-i-(i/ 2 )(2/) = *(~1)'<W^ H -l(i/2)(y) = -^YHf\ m {y). (55) 
However, for series of Bessel functions of the first and second kind, we have 

J-i-(im(y) = (-l)' +1 ^ +( i/ 2) (y), Y^_ {1/2) (y) = {-lfjf\ m (y), (56) 

that is, for v = 0, ±1, ±2, • • • the functions Ji+(i/2) and J_^_(i/ 2 ) (or, ^+(1/2) and Y_i_^/ 2 )) 
are linearly independent. In spite of this, by assuming that v 7^ 0, ±1,±2, ••• also for J 
and Y we guarantee that all of the solutions (160]) and ( 161]) for the spheroidal equation are 
two-sided since a* and 7* do not vanish for —00 < n < 00. 

On the other side, for two-sided series the domains of convergence are again given by ( 136]) 
with Ui substituted for Uj. For the one-sided series the solutions U j , in series of Bessel 
functions of the first kind, converge for any finite z. The behavior of the solutions at z = 00 
can be found from the fact that, for k fixed and |?/| — ^ 00 |17j . 



J «(y)~ y^cos^-f -£], Y K (y)~ ^±sm[y-f - f ] : | ar g2 /| < tt; 

H ^(v) ~ \ffv e^- f ] : 7r <C arg y <C 2tt; (57) 



H { l\y) ~ e-i y -^\ : -2vr < argy < vr. 
Now we consider the Meixner solutions. The substitutions 

y = l-2z, S(y) = z* [z-Vfiu{z) & S(y) oc [y 2 - 1]» U (z = ^f) (58a) 
transform the spheroidal wave equation (jHJ) into 

z(z - 1)^ + [- (// + 1) + (2/x + 2) z] ^ + [n (n + 1) - A + 4 7 2 z(z - 1)] [/ = 0, 
which is the CHE (pQ) with parameters 



^7T TV 



Zq = 1, B x = -//-!, 5 2 = 2/i + 2, 5 3 = /i(/i + 1) - A, w = ±2 7 , 77 = 0. (58b) 
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Instead of Z«(n), Meixner used the functions ip u J(v) which are given by [9] 

^(v) = ^^)zl% /2) (v), 

in analogy with the definitions of the spherical Bessel functions j K , y K , h . « and ft, 2 
by taking a; = —27 and using this notation, we get 



y + i 



J2 a ^nlMV ~ 1)], ^(j*, V) = s^c-^y), 



(59) 
. So, 

(60) 



o. 



1 _ 27(n+t/+l-/x)(n+i/+l) ol _ ( , \ / , , -■ \ _ > 1 _ 27(71+1/+//) (n+f) . 

™~ (2n+2v+3) ' 71 ^ )\ ) ' 'n — (2n+2v-l) ' 



and 



y+l 



El^r^tbdz + i)], CM = s? } (-/*,y). (61) 



(j), 



(2n+2f-l) 



?C3), 



For these solutions, conditions (154aj) and ( I54bl) reduces to 

2z/^0,±l,±2, •■• , z/± (//+ 1) ^ integer. 

The Meixner solutions are given by S [i ^{n,y) and S U 4 \fi,y). By the D'Alembert test S { j ] 
converge for |y — 1| > 2 (if % — 1,2) or for |y + 1| > 2 (if z = 3,4), as stated in (9j [19]. 
However, by the Raabe test they may converge at |y — 1| = 2 or \y + 1| =2 because relations 
(|36D and (pbjl yield 



>2, if 



Re(/i) < in 5?, 
Re(/i) > in S 1 ?; 



|y + 1| > 2 if 



Re (/i) < in 5 
Re (//) > in S 



U) 

3 > 

U) 



(62) 



(if Re (/i) = 0, the test is inconclusive). 

On the other side, the solutions w{u) for the Mathieu equation ([6]), considered as a 
particular case of the CHE, may be obtained by setting 



(63) 



w(u) — U(z), z = cos 2 (^), [a = l,i] =^ z = 1, 1 Mathieu eq. 
£i = -l/2, B 2 = \, B 3 = 2k 2 -&, ou = 4k, 77 = J as a CHE, 

where U(z) are the solutions for CHE with 77 = 0. Thus, from the previous solutions U^\z 
we get four sets of two-sided solutions w^(u) 

= J> ^ [4fccos 2 f], 



cos 



> 1, 



w% ] (u) = tan f £ i n + » + \) cos 2 f ] 



cos^f| > 1, 



(64a) 



w 



U) 



W = £(-1)X Z« [-4^m 2 f] 



; ( J (u) = cot f £ (-1)" (n + v + |) a n Z^ +1 h 4fc sin2 f ] 



sin ^| > 1, 



sin^ > 1, 



(64b) 



where the coefficients a n satisfy the relations 

2k(n + v + l)a n+ i+ a - 2k 2 - (n + + \) 2 a n + 2/c(n + u)a n -i = 0. (64c) 

For any solutions the behavior when z = cos 2 (au/ 2) — > 00 must be determined by using 
(1571) . If <t = 1 and u =real (Mathieu equation) the previous solutions are useless. Notice, 
however, that the one-sided solutions w ^ in series of Bessel functions of the first kind are 
convergent for all finite values of z. 
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2.2 Solutions for the Whittaker-Ince limit of the CHE 

To obtain solutions to the Whittaker-Ince limit of the CHE we apply [20] 

lim $ (a,c; -V) = r(c) 2/ ( 1 " c )/ 2 J c _ 1 (2 v ^), 

-me nc y^- c)l2 H^\{2^j), Im y > 0, (65) 



lim 

a— Kxi 



r(a + l-c) * (a,c;-^) 



on the hypergeometric functions used as basis for the expansions of the solutions for the 
CHE. For this it is necessary to rewrite the latter solutions in a suitable form and keep n 
fixed. Expansions in series of Y c -\ are obtained as a linear combination of the expansions 
in series of Hankel functions. In this manner, from the first set (119j) . we get a set of four 
solutions for Eq. (jlj). These are again denoted by U ^ (j = 1,...,4). In fact, we will 
compute only the limit of U\{z): the other solutions follow from the fact that the four Bessel 
functions satisfy the same differential and difference equations. 

On the other side, if U(z) = U(B\, B 2 , B 3 ; z ,q; z) represents an arbitrary solution for 
Eq. 01]), then other solutions are generated by the transformations ZT 2 and 37$ given by 

m 

3\U{z) = ^/^(d, C 2 , C 3 ; z , q; z), 

3%U{z) = (z- zof-^-^I^UiB^ D 2 , D 3 ; z 0} q; z), (66) 
3%U(z) = U(-B 1 - B 2 z ,B 2 ,B 3 - qz ;z ,-q; Zq-z), 

where C» and Di are defined in Eqs. ( 1T41 ). Thus, it is sufficient to take the limit of the 
first set of solutions ( f!9l) and of the coefficients (1201 . The other sets are obtained from U ^ 
through 

U® = 3T 2 U { [\ C/ ( ? = C/ ( f = W?- (67) 

First we find the four sets of solutions and use the Raabe test to study their convergence. 
In the second place, we write the solution for the Mathieu equation. 

2.2.1 The four sets of solutions 

To find the limit of U\(z) given in ( 1T91) . we rewrite that solution as 

qz 



Ef_1]i f \ az \n+v+l-(B 2 /2) 
1 ' ; m = $ 
V Om -4- 9;/ -4- 9 



n + v + 1 + irj, 2n + 2v + 2: 

17] 



(68a) 



T[2n + 2z/ + 2] 

where c n = b l n [-irf\- n and q = -2r/u. From a l n b l n+l + + 7n fo n-i = °> we S et 

- ir] alc n+1 + f3lc n + {-ir})~ l ^lc n -i = 0, (68b) 

where a\, fi\ and 7* are given in fT20]) . By supposing that n is fixed and using fl65]) . we obtain 
the solution U ^ in series of Bessel functions of first kind. In fact, we may verify directly 
that the four solutions given below satisfy Eq. (jS]). 
Then, the first set is given by 

V®{z) = z^l* 5>ir4zg +2i ,+i (2^qz) , (69a) 
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where, in the recurrence relations a!*c* +1 + /3*c* +1 + 7^c^_ x = 0, we have 

1 _ qZ0 [ n +v+2-^}[n+v+l-§±-^] 



a 



(2n+2u+2) (2n+2u+3) 



i~i \ / D , qz \B 2 -2]\B 2 + ^ 1 

# = B s -^+(n + i/+l-f) (n + v + f)- 2 J n+ J 2n+2 ;: 2) , 

! _ g^[n+«/-l+^][n+^+^ + ^] 



(69b) 
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(2n+2^-l)(2ra+2i;) 



The transformations (1661) . applied on L/y according to (IBTj) . generate the other sets, that is, 



^ - 



,i-b 2 - 



E(- 1 ) nc '^U + i (2v^) 



o; 



'i. 



(2n+2i/+2)(2n+2i/+3) 

^ = (,-, ) (1 ^ 2)/2 E(' - 

n 

g , [ w+ , +2 -^][ ra+t , +1+ ^ + ^ 
(2n+2i/+2)(2ra+2f+3) 



2n+2i>+ 
gzo[n+v+l--^-] [n+f 



7: 



1 ^| V 3 7 {j) 
L ) c n Zj 2n+2u+l 



T 3 
' /1 



2 *o 



(2n+2v-l)(2n+2v) ' 



2 - Zq) 



(2n+2v-l)(2n+2v) 



(70) 



(71) 



U { {\z) 



Z 2 o (z - Zb) 2 3 o 2 

n 



1 

v ) c n Zj 2n+2v+l 



2y^(z~^ 



z 



(/--() 



4_ ^o[n+,+l-^][n+,+^ + ^] 
7n (2n+2i/-l)(2n+2f) 



n (2n+2i/+2)(2n+2i/+3) ' 

For the these four sets of solutions, the conditions (1291) and (13"Uj) are replaced by 
2u, v ± ^ and 1/ ± + are not integers, 



(72) 



(73) 



while the relations ( 13"3"|) remain valid for the coefficients c^. 

The previous list completes the list given in Ref. [2j where the expansions U \ ' in series of 
J K and Y K have not been taken into account, whereas the expansions U^'^ have been written 
in terms of the modified Bessel functions K 2n+2u+ i [±2iy/qz\ and K 2n+2u+1 [±2i^ 'q(z — z Q ]. 
Moreover, now the regions of convergence are modified by the use of the Raabe test. 



2.2.2 Convergence of the solutions 

As in the case of the CHE, by the D'Alembert test the two-sided expansions U { { ] and U { f 



converge absolutely for \z\ > \zq\ while U % and U (J 4 ' converge for \z — Zq\ > \zq\. However 
by the Raabe test the solutions also converge at \z\ = \z \ and \z — z \ = \zq\ under the 
conditions similar (1361). that is, 



(./) 



\z\ > \z \ if 



fRe 


B 2 + 




< 1 in 






20 




1 Re 


B 2 + 


Bi 


> 1 in 




ZQ _ 





rU) 



r(i) 
2 5 



\z — Zq\ > \zq\ if 



and 


\ z 


~ Zq\ = 


fRe 


Bt 


> -1 




. z ° . 




[Re 


'Sx 


< -1 




. z ° . 





rU) 

3 ) 

r(i) 

4 ' 



(74) 



The test does not assure convergence at z = 00 and, so, the behavior at z = 00 again 
deserves special attention. We will find that the one-sided infinite series U ^ in series of 
Bessel functions of first kind converge for any finite value of z. 
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Relations flH} correspond to (136)) with the replacements -H- Uj. In fact, a few 
modifications in the previous demonstration lead to (1741) . Thus, if n — >■ ±00, we find 



Si 

2(1 



- (b 2 

When n — > 00 the minimal solution for c n+ \/c n is 

1 - i f2l/ - 5 2 



£- + [4n(n + 2i/ + l)] + 
^ = 0. 



5i 

«0 



4n 2 

qzo 



2v - Bo 



Mi 

zo 



(75) 



while the minimal solution for c n _i/c n , when n — > — 00, is 

2v + B 2 



r 1 

"n+1 



An* 
qzo 



§1 

ZQ 



(76) 



On the other side, the behaviors (1A.21j) and (lA.22j) for the Bessel functions lead to 



lim ^ 2n+2u+3 

n^oo J 2n+ 2v+l 



lim 



lim 



7(2,3,4) 

; 2n+2^+3 _ 

7(2,3,4) 

y 2n+2i/+l 

7 (j) 



qz 
An 2 

An 2 
qz 

An 2 



1 + 



1 

n 

1 



2z/ + 



2v 



Thus, when n tends to +00 



qz 



1 

n 



2u 



[j = 1,2,3,4]. 



(77) 



lim 

n— ¥00 



lim 

n— >oo 



Cn+ 1^2n+2iH-3 


g 2 z z 


~l 7O) 
Si+l Zy 2n+2i/+3 


1672 4 
= -*> [l 

z 


G n Zy 2n+2i/+l 



- ( 4z/ + 6 

n 



1 - 



a-* 



Bo 



Bi 

z 



[J = 2,3,4] 



(78) 



and, when n — > — 00, 



lim 



Cj) 

2n+2v-l 



L n Zy 2n+2i/+l 











1 H 


-i( 


z 




n \ 



\-B 2 -^ 

z 



[j = 1,2,3,4]. 



(79) 



Hence, by the Raabe test the expansions are convergent for \z\ > \zq\ as indicated in (1741) . 
For the other sets of solutions the domains of convergence follow from the transformations 
(16"6"j) applied on U ^ according to ( 167ft . Moreover, from the first limit given in (I78p we see 
that the one-sided infinite series U j converge for all z excepting probably the point z = 00. 
The behavior any solution when z — > 00 must be determined by using (1571) . 
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2.2.3 Heine's solutions for the Mathieu equation 

From the previous solutions we recover the usual solutions in series of Bessel functions 
for Mathieu's equation (called Heine's solutions [21]) by means of the substitutions 

w(u) = U(z), z = cos 2 (au) z = 1, 1 Mathieu eq. as Whittaker- 

B x = -\, 52 = 1, £ 3 = £-£, q = k 2 J Ince limit of the CHE. 

However, the regions of convergence for some solutions turn out to be improved by the Raabe 
test. This fact is useful for some applications, as we will see in Sec. 4. 
Relations fl33|) . with b l n replaced by d n , yield 



cl = {n + v + |) 4, c\ = {-l) n cl 4 = (-IT {n + u + \) c\. 
Then, by writing c n = c\, w^(u) = U^ (z) and setting y/k 2 = k we find 

w V>(u) = ("l) n Cn4nU+i [2kcos(au)} , I C0S(<7U)| > 1, 

n 

+{\u) = tanH ^ (-1)" ( n + i/ + |) c n zg +2 , +1 [2fc cos(a M )] , | cos(au)\ > 1, 



f81a) 



w 



- c n Z^ +2u+l [2ki sm(au)\ , | sin(a M ) | > 1 

= cot((7u) (n + v + |) c n Z 2 ( J f2i/+1 [2ki sin(an)] , | sin(o"w)| > 1, 



(81b) 



where the coefficients c n satisfy the relations 

k 2 c n+1 + [(2n + 2u + l) 2 - a] c n + k 2 c n ^ = 0. (81c) 

As in the case of the two-sided solutions f)64ap and (164bj) . obtained from the CHE, the above 
solutions are useless for a = 1 and u =real (Mathieu equation). 

The conditions ( !73|) reduce to 2v ^ Z and assure the linear independence of the terms in a 
given series. In addition, the above notation for solutions of the Mathieu equation is similar 
to the one used by Erdelyi [21]. However, in Ref. [T7J Q2] the coefficients (c n +i, c n , c n _i) 
are replaced by (c 2n+ 2, c 2n , c 2n _2) and the Bessel functions Z 2n +2v+i written as Z 2n +v'- this 
is equivalent to put 2v + 1 = v with v ^ Z. The above domains of convergence may be 
compared with the ones of solutions 28.23.2-28.23.5 of |18j . 

By the Raabe test, the the two-sided solutions w^{u) and w^{u) are absolutely con- 
vergent also at | cos(<7w)| = 1 and | sin(<TM)| = 1, respectively. In Refs. [7J [T7J [191 122] these 
points are not included in the domains of convergence due to the use of the D'Alembert 
test. The one-sided solutions converge for any finite u; in Ref. [IB] it is stated that this 
property holds also for two-sided solutions (a misprint, for certain). 



2.3 Solutions for the double-confluent Heun equation 

For zq — > the limits of the solutions of the sets ( JT9"j) and (|2"3"|) are equivalent, respectively, 
to the limits of (125]) and ([27]). This gives a subgroup constituted by two sets of solutions. 
A second subgroup, with different regions of convergence, results from the transformations 
of the DCHE and cannot be obtained as limits of the expansions in series of Coulomb wave 
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functions. In effect, if U(z) = U {B\, B2, B 3 ; 00, rj; z) denotes one solution for the DCHE, new 
solutions can be generated by using the transformations [3] 

hU(z) = e^z^^Ui-B^A -B 2 ,B 3 + 2- B 2 ; uj, 77; z), 

t 2 U{z) = e^+^z-^-^U (B 1 ,B 2 ,B 3 ;u = 1, 77; f^) , (82a) 

t 3 U(z) = U(B 1 ,B 2 ,B 3 ;-u,-r l ;z), 

where, on the right-hand side of the second relation, we have 

B 1 = uB it B 2 = 2 + 2it], B 3 = B 3 - [f + ir)] [f - ir) - l] , if) = f - 1. (82b) 

The transformation t 2 is the one responsible for the second subgroup of solutions. 

In the following we use the notation given in f[T7|) and ffTS]) also for solutions of the DCHE. 
We find that the Raabe test is useless for the DCHE. 

2.3.1 The two subgroups of solutions 

The limit zq — > does not change the form of (fl~9j) but modifies the series coefficients. 
Thus we have 

b 1 \2ituz\ n+u+1 ~ : ^~ 

U x {z) = Yl no 4.9 T 01 $ {n + v + l + iv,2n + 2v + 2; -2iuz) , 



T[2n + 2u + 2] 
b 1 \—2i(jJz\ n+y+1 ~^2' 

Uf(z) = e ±iujz V i (n + u+l± ir), 2n + 2v + 2; T^iuz) 

T[n + v + 1 =f ir]\ 



f83a) 



whereas coefficients (|20|) when zq — > give 

l 2iujB! [n+i/+2-^-] 

a n ~ (2n+2i/+2)(2n+2i/+3) ' 

ft = B 3 + (n + v + 1 - f ) (n + v + f ) + (2ra 2 ;gg^ 2) , (83b) 



^ 2iuiB\ |ra+i/+-^ — lj (n+^+j»7)(n+i/— ir;) 
(2n+2i/-l)(2n+2i/) 



On the other side, since 



lim f 2 (z) = lim 

ZQ-»0 2 — >0 



the solutions obtained from (T53D and (1241) are 



n+iH-1 ^ 



(84a) 



t/- 2 ( z ) = e iuz+ ^ V ^ — $ fn + v + 1 + Z77, 2n + 2zv + 2; -2iuz] , 

v ; ^ r\2n + 2v + 2] 1 1 > J> 

s A 2 [_ 2iu z] n+u+1 ~ 

U±(z) = e ±iuJZ+ ^ > - -tf [n + v + 1 ± ir), 2n + 2v + 2; =F2iwz] , 

r n + 1/ + 1 =F iri\ 

n L ' 

where the coefficients for the recurrence relations are given by [ft — ft) 

2 —2iu)Bi[n+u+^~^ 2 2iu)Bi\n+v+l—^]{n+i/+iri){n+v—iri) (RAh\ 

a n ~ (2n+2i/+2)(2n+2iy+3) ' _ (2n+2i/-l)(2n+2v) ' ^ ' 
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Notice that U2 = t\Ui. We can check that £3 does not generate new independent solutions 
when applied on Ui and U2. Furthermore, by using i 2 we construct a second subgroup having 
two sets of solutions. Thus, we form a group of four sets of solutions Uj, given by 

UiC*), U 2 (*) = hU^z); U 3 (z) = hViiz), U 4 (s) = t 2 U 2 (z) = t 3 V 3 (z) (85) 

where U; = (U h U+, U~) as in Eq. OTjl. 
The solutions of the third set are 



UaW 2^ r[2n + 2i/ + 2] 



-Bo B\ 
n + v + —,2n + 2v + 2; — 

2 z 



\r /■y] n+u+ ^- (86a) 

* [n + 1/ + 1 ± 2n + 2u + 2; ±^] , 
where, in the recurrence relations for (/3^ = 

3 _ 2jt£gi (n+i/+l-iT)) 3 _ 2«jBi [n+t/-l+^] [n+y+1-^-] [n+u+irj] 

a n ~ (2n+2u+2)(2n+2u+3) ' ~ (2n+2i/-l)(2n+2i/) 



(86b) 



Finally, since the fourth set is given by 

[CM*), = faEM*),**^)] = [t 3 U 3 (z),t 3 Ut(z)] , (87) 

it is obtained by replacing (77, u) by (—77, —a;) in (U 3 , U^ 1 ) and respective recurrence relations. 
For these solutions of the DCHE, the restrictions (129]) are replaced by 

2v and v ^pin cannot be integers in Ui and Uo , 

(88) 

2v and v =p cannot be integers in U3 and U4 . 

Under these conditions, the three hypergeometric functions of a given set are linearly inde- 
pendent and, consequently, each solution can be written as a combination of the others by 
means of (1A.5|) . However, to have two-sided infinite series, in addition (IHHj) . we must require 
that 

v =1= ^ cannot be integer in Ui and U2 , v =F 277 cannot be integer in U3 and U4 (89) 

These conditions take the place of ( 130]) . Eqs. ( ]3Tj) hold also for the present solutions, but 
Eqs. ( 133]) are replaced by 

_ [_i]»r[n+i/+2-^] ! 3 _ r[n+E/+2-g&] x 4 _ [-l]"r[n+i/+2-^] x 

°«- r [n+i/+^-] n ' n r[n+u+l-i V ] °ni °n r[n+u+l+ir)} °n- l JU J 

2.3.2 Convergence, Whittaker-Hill and Mathieu equations 

The two subgroups present different regions of convergence: 

• The solutions Ui and U2 converge absolutely for \z\ > but may be unbounded at 
z = 00. The one-sided series (n > ), U\ and U 2 , converge also at z — 0. 

• The solutions U3 and U4 converge for \z\ < 00 but may be unbounded at z — 0. The 
one-sided series, U 3 and {7 4 , converge also at I2I = 00. 
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Once more we take the solutions ( I83al) in the form ( jl~2l) . Then, the ratios 6* +1 /6* and 
b\_ x jb x n are obtained from 



n 



bl 



2n 



n + 2v + 1 



iuBi'n 



bl 



which takes the place of Eq. (1381) . The minimal solutions are 

h\l_ 2^+2 1 . b} 

n L n. J 



rt — > oo : 
n — > — oo : 



fe n+l ^ C 



6J, 



2in 



n-i 2(n+2i/+l) 



„+l 2n 2 (n+2;/+l) 
— 1 — rvj 1 - 



iujBi 



Combining these limits with (j40p . we find the relations 



lim 

n— >oo 



^ +1 n+I ,+l(^,^) 










4ra 3 



lim 



K^n+viV^z) 





B 1 




zn 



(91a) 



lim 



b x n+1 ^ n+u+1 (r],ujz) 



B 1 



zn 



(91b) 



which follow as well from (141 ap and (141bp as z — > 0. Thence, by the ratio test, the solutions 
Ui converge for \z\ > \Bi/n\, that is, for \z\ > 0. The same holds for the second set of 
solutions, U2 = tiUi. For 2—7-00 the solutions may be unbounded insofar as their behaviors 
are formally given by f|44|) and (j4"5l) . The corresponding one-sided solutions in series of 



regular confluent hypergeometric functions converges for \z\ > because the limits (191bj) 
become irrelevant. 

On the other side, by applying t 2 on ( 191 aft and ( I91bl) we find 



lim 

n— ¥00 



lim 



n Yn+u 



ujB\ 



in 3 z 



lim 

n— ¥00 



6» + , ^) 



2ujz 



n 



2ujz 



n 



(92a) 



(92b) 



Then, the solutions U3 converge for finite values of z\ the same is true for U3 (if n > 0, 
U3 and U4 converge at z = 00 as well). However, we must determine the behavior of each 
solution at z — 0. Thus, from relations (IA.7P 



lim ^3(2) 

Z->0 



E 



6 3 



1> + v + 2 - B 2 /2] 



z 

By 



2-B 2 



r_l\n+u+B 2 /2 £3 

T[n + u + B 2 /2] 



(93) 



while from Eq. (I A . 6 1) 

E 



lim (7+ z 

z— >0 



lim C/q 



r[n- 



62 



5 a /2] 

1/2 E 



-1 



+B2/2 £3 



T[n 



B 2 /2\ 



3tt 
2 



3tt 
2 



< arg (^) < 

< arg (-^l) < 



3tt 
2 ' 



3tt 
2 • 



(94) 



Therefore, U3 and t/ 3 may be unbounded at z = as occur, for example, if Ke(B±/ z) — > 00. 
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Now we consider the properties of the solutions for the Whittaker-Hill and Mathieu 
equations obtained from the solutions of the DCHE. The Whittaker-Hill equation ([7J (WHE) 
is transformed into the DCHE (j2J) by the substitutions 

z = e 2 -, W(u) = z^^e^U(z) =► B ± = - J, j WHE as 

B 2 =p + 3, #3 = (i + l) 2 + fi-f, iw = f, <i7 = j(p + l) J a DCHE. (95) 

So, by inserting the above expressions into the solutions Ui(z) = {Ui{z) 1 Uf{z)) for the 
DCHE, we get the solutions W;(w) = (Wi(u),W^(u)) for the Whittaker-Hill equation. Thus 



Wl = e , tP ,«,.. /+ £ [„ + „ + EH, 2n + 2 „ + 2; _£ e -»] , 

n 



(96) 



where 
and 



y+ = f + (u) = e ^l A ) cos ^ u \ f- = f~(u) = e*^/ 4 ) sin ( 2?u ) 5 



w} = e ^ +ll v £ *j^gj> [„ + „ + 1 ± £±1, 2„ + 2, + 2; T $«*■] 



(97) 



where 



g+ = g+( U ) = e i(e/4)sin(2,„) ) ^- = ^- (u) = e _ K /4) cos(2 fM ) 

The other sets are given by 

W 3 (u) = Wi(-u), W 4 («) = W 2 (-u). (98) 

In the recurrence relations ( 1T51) for and 6 2 , the coefficients are obtained by replacing the 
parameters (1951) into (I83bj) and (I84b|) respectively. The convergence of the two-sided infinite 
series is given by 

\e 2i " u \ > for Wi(«) and W 2 (u), \e~ 2iiU \ < oo for W 3 (u) and W 4 (u). (99) 

Then, if <j = % and u is real (WHE, not modified WHE) the four sets of solutions are 
convergent for all values of z = exp (2iu), in opposition to the solutions obtained from the 
solutions for the the CHE - see Eqs (1571) . 

On the other side, the Mathieu equation (jfJJ) is converted into a DCHE by the substitutions 

z = e lau , w(u) = z 1 l 2 e ik l z U(z) ] Mathieu equation 

B t = -2ik, B 2 = 2, £ 3 = (1/4) -a, u — k, rj = J as a DCHE. 

For these parameters the preceding solutions Uj(z) for the DCHE are expressed in series of 
Bessel functions J re , .Hy and « by using the relation (|A. 19[) . We can as well include the 
function Y K since these are linear combinations of the Hankel functions. Thus, from Uj(z) 
we obtain the sets of solutions, W;(u), 



(101a) 
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where the recurrence relations are 

2k " [SfSd «n + i +[\-*+{n + u)(n + v+l)] a n + 2fc 2 a n _ x = 0. (101b) 

The convergence of the above series requires that 

\e iau \ > for w { {\u) and w%\u), \e~ iau \ < oo for w { {\u) and w { $(u). (101c) 

2.4 Solutions for the Whittaker-Ince limit of the DCHE 

For the limit ([5]) of the DCHE we find only three sets of solutions. Two of them are 
obtained by applying the Leaver limit (zq — > 0) on the solutions in series of Bessel functions 
( I69aj) and (1701 for the Whitakker-Ince limit of the CHE; the other set is given by series of 
confluent hypergeometric functions and is obtained by applying the Whittaker-Ince limit on 
the solutions fl86al) or flHTj) for the DCHE. We have found no transformation of variables 
connecting the solutions in series of Bessel functions with the solutions in series of confluent 
hypergeometric functions. 



2.4.1 The two types of solutions 

When z — > the set (I69al) retains its form, that is, 

U { {\z) = z^ B ^ ^-^"c^SU+i (102a) 
?! 

but the coefficients ( I69bl) for the recurrence relations reduce to 

n l = gBi[n+*+2-ga] 
"n (2n+2i/+2)(2n+2y+3) ' 

= B 3 + (n + »/ + 1 - f ) (n + i/ + f ) - (2n+ f;)(S;£ +2) , (102b) 

1 _ g Bi[n+v+^-l] 
~~ (2n+2u-l)(2n+2u) ' 

This set arises also from ( 17TT) when Zq 0. It can as well be obtained by taking the 
Whittaker-Ince limit of solutions (183aj) provided we add the solution in series of Y2 n +2v+i- 
On the other hand, when zq — > the sets ( ITU]) and ( !72|) give 



= e^/V 1 -*)/" ^(-l)"4zS +2 , + i (2v^) , (103a) 

with 

£2 = ol 2 = _ qB^n+u+l-^f.) , 1Q3b) 

rni In (2n+2v-l)(2n+2u) ' v y 



o: 2 



ij_Bi(n+i/+-^) a2 _ ol 2 qBi(n+v+l- 

n ~ (2n+2i/+2)(2n+2iH-3) ' ~~ ^n) 7n — — (2n+2u-l)(2n+2u) ' 



(2) 

in the recurrence relations for c n ■ This set can also be obtained by taking the Whittaker- 
Ince limit of solutions f !84ap . On the other side, solutions (j86aj) and 087|) yield the same 
Whittaker-Ince limits in series of confluent hypergeometric functions, namely, 

= £ 'rf^Cf t [» + , + f , 2» + 2, + 2; f] , 

(104a) 

M = £ rg^gr.^ * [» + " + 1 ± V. 2n + 2, + 2; ±&] , 



?[n+i/+lT(B2-2)/2] 
n 
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where, in the recurrence relations for , 

-qBi 



(2n+2u+2)(2n+2u+3) ' Hni In (2n+2u~l){2n+2u) " v ' 



This set yields expansions in series of Bessel functions if B 2 = 2. 
The three sets of solutions are given by two-sided infinite series if 

2v and v ± ^ are not integers. 



(105) 



The condition 2^ ^integer avoids vanishing denominators in the recurrence relations and 
assures the independence of the Bessel and hypergeometric functions in each series. Besides 
this, if v ± (B 2 /2) is not integer, the coefficients are connected by 



■1) 



n r[n+y+2-(B 2 /2)] l 
T[n+u+{B 2 /2)\ C n' 



C 3 = r[n + !/ + 2-f 



(106) 



For completeness we mention that, if U(z) = U(Bi, B 2 , B 3 ; z) is an arbitrary solution 
for the Whittaker-Ince limit of the DCHE, then a new solution may be generated by the 
transformation T whose effect on U (z) is 



TU(z) = e B ^ z z 2 - B2 U{-B x A -B 2 ,B 3 + 2- B 2 ; z) 



(107) 



We can check that TU^ = t/j while T(U 3 , Uf) does not produce new solutions. 

The solutions in series of the Bessel J 2n + 2v+ \ and Y 2n+2u+ i, as well as the solution U 3 (z) 
in series of regular confluent hypergeometric functions, are new, whereas the other solutions 
have already appeared in [2]. 

2.4.2 Convergence of the solutions 

As in the case of the solutions for the DCHE, we find that 

• The solutions and U^ 2 converge absolutely for \z\ > but may be unbounded at 
z = 00. The one-sided series, U ± and U21 converge also at \z\ = 0. 

• The solutions (Us, U 3 ) converge for \z\ < 00 but may be unbounded at z = 0. The 
one-sided series U3 converges also at \z\ = 00. 

In effect, if n — > ±00 the recurrence relations for cL give 



qBl c n+i 



Then, the minimal solutions are 



[4ra(ra + 2z/+l)] - afi^i = 



n — > 00 : 
n — > —00 : 



C l 4(n+l) 2 {n+2u+2) ~ An* \_ 1 ^ W \n)i 

c n-i ^ qBi _ qBi , n ( 1\1 

c\ i{n-l) 2 {n+2v) ~ 1^ I 1 "r W \n ) J 



Thus, by using the ratios ( I77|) we find 



lim 

n— ¥00 



lim 

n— >— oo 



-\J2r1 



+2i>+3 



clJ: 



2n+2u+l 



q 2 B lZ 
16n 5 



2n+2v-l 



^n Zj 2n+2v+l 



Bi 

nz 



lim 

n— >oo 



r- 1 

L n+1 



„1 7 (2,3,4) 
c n+l Z/ 2ra+2t/+3 
„1 7(2,3,4) 
G n Zy 2n+2;/+l 



4n 2 (n+2z/+l) 

4ra 2 (n+2t/+l) 
9B1 



-5! 



— , J = 1,2,3,4 



108) 



(109) 
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[these expressions also result as limits of (ITS]) and ( 1791) when z 0]. Therefore, by the 
D'Alembert test the two-sided series solutions converge absolutely for \z\ > 0; for one- 
sided series, U \ converges also at \z\ = because only the limit n — > +oo is relevant. 
However, the behavior of each solution at z = oo must be analyzed by using the relations 
(I57j) . The convergence of the set is inferred from = TU since the transformation 



T replaces Bi by —B\, it does not change the convergence properties for the series in 
but can change the behaviors at z = oo - see Eq. (I107j) . 

To get the convergence of the third set, we write the solutions ( I104ap as 



U 3 (z) = z ^e B ^ 



i(2-B 2 ) iB 1 



2 7 2z 

where we have used the definitions (IA.10I) for ^ n+u . When n — > ±oo, (? n obeys 

a* £|l _ [4n(n + 2v + 1)] - qB^ = 
whose minimal solutions are 

c n+l -qBi qBi \-i , n ( 1 \ 1 4-1 4n(n+2i/+l) r . 1 
-^T 4(n+l)(n+2„+2) = I 1 + \n)\ ~ n Wl ' V 1 ~* °°\ 

c n-l qBl qBi N i n I \\\ _ 4+1 4n 3 (n+2v+l) r 1 

~ 4(n-l)*(n+2,/) = 4^ I 1 + \n>\ & ~| > [7l->-OOj. 

On the other hand, for finite values of y = iB\/{2z) relations ( 1A.16I) and ( 1A.17I) give 



(110) 



4>n+v 



4n 2 z 



Ill) 



Thence 



4 + l«W„ + l 


~ 9 l 


C n + l^n + i/+l 




[n —7- oo] 


C 3 0„ +I/ 


16n 4 z' 






'n—l^rt+f—l 




^n+u = ( 0n- 


7/; (±) ^ 
hi') rn+i/ 1 


[n — > — oo 



;ii2) 



Thus, the two-sided series solutions U 3 (z) converges for finite values of z, excepting possibly 
the point z = 0. Since the limit n — > oo is the only one applicable for one-sided series, the 
solution U 3 converge also at z = oo. In order to verify if the solutions are bounded or not 
z = 0, we must find the behavior of such solutions at z = 0. The results are the same given 
in Eqs. ([93]) and (El for the DCHE. 



3 One-sided Series Expansions 

In Ref. [TU] one-sided solutions in series of Coulomb wave functions have been obtained 
from two-sided solutions. In this section this problem is reconsidered because the redefi- 
nitions of the Coulomb functions modify the recurrence relations for the series coefficients, 
while the Raabe test modifies the regions of convergence for the solutions of the CHE. We 
examine also the Whittaker-Ince limits of the CHE and DCHE and the conditions for having 
finite-series solutions. 

In fact, it is possible to restrict the summation of the two-sided series to n > by writing 
the parameter v as function of the parameters of the CHE. For each of the four sets Ui(z) 
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we get two expressions for v and, accordingly, eight sets Uj(z) of one-sided series solutions. 
Thus, for UiO) we obtain v = (B 2 /2) - 1 and v = (B 2 /2) + (B x /zq) - see Appendix B. We 
consider only the set Ui(z) obtained by setting v = (B 2 /2) — 1 in Ui(z), and generate the 
other sets by applying the transformations ({TBI) in the order given in Eqs. (11141) . If rj ^ 0, 
there are three possible types of recurrence relations for the series coefficients; if r\ = 0, only 
two types. 

Solutions for the DCHE, as well as for the Whittaker-Ince limit of the CHE and DCHE, 
are obtained by the procedures used in section II. We will see that in a fixed set for the CHE 
and DCHE, it is possible to have one solution given by infinite and one given by finite series. 
On the other side, the one-sided expansions in series of regular hypergeometric functions 
<3>(a, c; y), their particular and limiting cases, are convergent also in the vicinity of z = 0, as 
we have already seen in section II. 

First we regard the solutions for the CHE and DCHE; after that, the solutions for the 
Whittaker-Ince limit for the CHE and DCHE. 

3.1 Solutions for the confluent Heun equation 

The one-sided solutions Ui(z), obtained in Appendix B, are 

V 1 {z) = z-%j£blV n _ 1+ ^{r,,u,z), = , (113) 

n=0 ^ 2 ' 

which can be rewritten as in Eqs. fll 22j) . Then, the transformations ([TBI give the sets: 

Ui = Ui|_ , U 2 = TiUi, U 3 =T 2 U 2 , U 4 = TiU 3 : first subgroup; 
U 5 = T 4 Ui, U 6 = T1U5, U 7 = T 2 U 6 , U 8 = TiUr : sec. subgroup. 

We can check that these solutions are obtained by setting v = vi in the two-sided solutions 
Uj (and also in a l n , /?* and 7*) according to 



Ui=UJ , U 2 = U! , U 3 = U 2 , U 4 = U 2 ; 

u 5 = u 3 | , u 6 = u 4 | , u 7 = u 4 | , u 8 = u 3 | 



(115) 



with 



^l = f-l, ^2 = f + f, ^3 = l"f, ^4 = -f-f, U iH = U it (116) 



where i — 1,2, 3, 4. Thus, v is given by the same expression in the solutions Uj and Uj +4 . 

First we write the general form of the recurrence relations for the series coefficients and 
state the convergence of the infinite series. Then, we examine the conditions for having finite 
series for the cases 77 7^ and 77 = 0. Actually, here we write only the sets Ui and U5; the 
other sets are given in Appendix C. 

3.1.1 Recurrence relations, transformations and convergence 

In the following a n , (3 n , % and b n mean ct l n , (3 l n , 7^ and b l n . Then, by Appendix B the 
forms of the recurrence and characteristic relations are 

otobi + /3 b = 0, I o « o7l a l72 a 2 j 3 

> =>• Po — — 5 3 — • • • , ( 1 1 / ) 

a n b n +l + Pnb n + 7n&n-l = 0, U > 1 \ ft- ft- ft- 
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«o^i + $obo = 0, aib 2 + 
Pik + [71 - {rf + |) d-j] 6 = 0, 
a n b n+1 + f3 n b n + 7„6 n _i =0, n > 2 



/So 



«o [71 - {v 2 + |) «i 



72 



'1181 



$0 — irja- 



bo = 0, 



a n b n+ i + f3 n b n + jnbn-! = 0, n > 1 



/So — ^eL 



«07l «l72 «273 

4i- A- &- 



(119) 



where a n , $ n and 7 n depend on the parameters of the equation. Notice that n > — 1 for a n , 
n > for /3 n , and 77 > 1 for 7 n . 

Besides this, relations (I117H119P stand if the series really begin at 77 = 0, that is, if 60 7^ 0. 
If the summation begins at 77 = N > 1, the recurrence relations are given by (1128aj) . For 
series beginning at n — 0, the recurrence relations for 77 7^ and 77 = are given by 



rj^O: ((EE} if 2ut + 2 ^ 1,2; (TTTHD if 2^ + 2 = 1; (THID if 2z/ t + 2 = 2, 
7/ = : (mZD if 2i/i + 2 ^ 1; (EH]) if 2^ + 2 = 1 



(120) 



where Vi are the parameters given in f )116p . Thus, U, and Uj+4 {% = 1, • • • , 4) have the same 
type of recurrence relations. On the other side, to assure the linear independence of each 
term of the series we demand that 



2b>i + 2 7^ 0, — 1, —2, • • • for series beginning at 77 = 0. 



'1211 



These restrictions also avoid vanishing denominators in the recurrence relations. Explicitly 
the initial set reads 



Ui(z) 



U?(z) 



^ bl [2iuz] n ^ ( B 2 

^ „ ( 77 + it] + 2ti + B 2 ; -2iuz 



D ±iuiz 



bl \-2iuz\ n 



£ 

n=0 



Bo 



* n ± irj + — , 2n + 5 2 ; +27wz , 
r[77 =f 777 + {B 2 /2)\ \ 2 



(122) 



where the recurrence relations for b\ are given by (11261) with 



2iuizq [77 + 1] 



a, 



7? 



(2n + B 2 ) (2n + B 2 + V 



uvzo [-82+^1 [B2-2] 

$ = - n ( n + B 2 - 1) - wt/zq - B 3 - (2n ^ 2 _ 2 )° (2 J w+B2) : 



2iuz [n + B 2 - 2] 



7« 



77 + 52 + ^-1 



2() 



[n - 1 + 77? + ] [n - 1 - «77 + ^] 



(277 + B 2 - 3) (277 + £ 2 - 2) 



(123) 



The fifth set, U 5 (z) = T 4 Ui(z), is 



^ [27a;(z - z )F ^ 



n=0 



T[2n + B 2 ] 



B 2 

n + ir]+ — , 277 + £? 2 ; -2mj(z - z . 



n=0 



r[77 T 77/+( J B 2 /2)] 



B 2 , . 

77 ± 777 + — , 277 + B 2 ; ^2iu{z - z ) 



(124) 
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with 

2iuizq [n + 1] 



7 5 



n + B 2 + ^ 



6 5 = 1 



(2n + B 2 ) (2n + B 2 + 1) 
2iuz [n + B 2 - 2} n - 1 - [n - 1 + irj + ^] [n - 1 - irj + 



(125) 



(126) 



(2n + B 2 -3) (2n + 5 2 - 2) 

The recurrence relations for 6* and are given by 

rj ^ : Eqs. (EZD if 5 2 ^ 1, 2; Eqs. $111) if S 2 = 1; Eqs. (EH]) if B 2 = 2. 
r/ = 0: Eqs. (HIZD if 5 2 ^ 1; Eqs. (HH if B 2 = 1, 

while for both sets conditions ( 1121ft become 

£ 2 ^ 0, -1, -2, ■ • • for Ua and U 5 . (127) 

Although the recurrence relations f II 1 8 j) and f II 1 9 j) hold only for special cases, their rele- 
vance becomes apparent by regarding the Whittaker-Hill equation as a CHE. By Eqs. (1511) . 
for this case Zq = 1, Bx = — 1/2 and B 2 = 1. Then, examining also the solutions written in 
Appendix C, we find: Eqs. (TTT8D for Ui and U 5 ; Eqs. (TTT7j) for U 3 and U 7 ; Eqs. (1TT9|) for 
the other solutions. 

The relations (I117H119]) and (11211) hold if the series begin at n = and extend to infinite. 
However, in some cases the series truncate on the right-hand side leading to finite series (see 
page 146 of [7j): for these solutions the ratio tests are meaningless. In other cases the series 
are truncated on the left-hand side and, then, the series begin at n > (page 171 of [7], Ex. 
2). 

In effect, if 7„ = for some n — N > 1, we have three-term recurrence relations only if 
the series terminate at n — N— 1, that is, if b N = b N+ i = ■ ■ ■ = 0. In particular, for N = 1 the 
series presents just the first term except when the recurrence relations are given by ( 11181) : 
in this case the series presents only the first term if — [rf + (l/4)]a_! = 0. 

On the other hand, if a n = for some n = N — 1 > 0, then the theory of three-term 
recurrence relations remains consistent only for series beginning at n — N, that is, if bo = 
■ ■ ■ = 6 N _i = 0. In this case, the recurrence relations and the corresponding characteristic 
equation (in terms of continued fractions) are 

a N fr N +i + /3 N fr N = 0, o d N 7 N+ i d N+ i7 N+2 , , 

> => Pn = ;= • (128a) 

a n b n+1 + f3 n b n + %b n -i = 0(n>N+l) J p N +i~ Pn+2~ 

In this case, conditions f 1 1 2 1 j) are replaced by 

2ui + 2N + 2 ^ 0, -1, -2, • • ■ for series beginning at n = N. (128b) 

The convergence of the one-sided infinite series is obtained from the convergence of the 
two-sided infinite series. We find that 

• The expansions Ui(z) in series of regular hypergeometric functions converge for finite 
values of z in the complex plane, as explained in sec. II. 

• By the D'Alembert test the expansions Uf{z) in series of irregular hypergeometric 
functions converge for finite values of z such that \z\ > \z \ (if i — 1, • • • , 4) or \z — Zq\ > 
\zo\ (if i = 5, • • • ,8); by the Raabe test, they converge also at \z\ = \zo\ and \z— z \ = \zo\ 
if the conditions given in (11291) are fulfilled. 
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• The behaviour of the solutions at z = oo must be determined by using ( 1A.6I) and ( 1A.7fl . 

The regions of convergence given by the Raabe test result from the combination of relations 
(EED with (flT6l) . So, we find 

> -i = u±; 

(129) 

< "I = U± 7 , 

where U^ n means: and U^. The test is inconclusive if Re[£>2 + (Bi/zo)] = 1 and 

~Re[Bi/z ] = — 1 in first and second cases, respectively The Raabe test is not valid if the 
above restrictions are incompatible with the restrictions which assure the independence of 
the hypergeomertic functions in each series; in addition, it does not assure that the solutions 
are bounded at the singular points. 



\z\ > \z \ if 




<1: U±; 



> 1 



Zo\ > \ z o\ if 



|Re 


Mi 






1 Re 


Bi 




ZQ 



3.1.2 Finite-series solutions for rj ^ 

If rj 7^ 0, in the same set it is possible to have solutions in which the series truncate on the 
left-hand side (this gives infinite series solutions) and on the right-hand side (finite series). 
In this case, the finite and infinite series present different series coefficients, despite the fact 
that these are formally denoted by the same symbol b l n . 

We consider only the first set but the following considerations can be extended for the 
other sets. Thus, if 

B 2 

2n + B 2 and n ± in + — are not zero or negative integers in Ui, (130) 

for any admissible n, each of the solutions (11221) can be expressed as a linear combination of 
the others by means of Eq. flA.5j) . For series beginning at n — 0, the above conditions are 



satisfied if B 2 , irj + B 2 /2 and — ir] + B 2 /2 are not zero or negative integers. However, when 
the second or the third condition is not valid, the solutions can be interpreted as (rj ^ 0): 

• If ir) + B 2 /2 = —£, where i is zero or a natural number, then Ux and (jf are solutions 
linearly dependent given by a finite series with < n < i because 7^ = ^ +1 = 0, whereas 
the series in fTf begins at n = i + 1 because \jY{n — £) — for n < t: in this case the 
solutions cannot be related by (IA.5j) . 



If —irj + B 2 /2 = —£, where i is zero or a natural number, then JJ\ oc U 1 is given by a 
finite series with < n < £, whereas the series in begins at n — £ + 1. 



The linear dependence follows from the generalized Laguerre polynomials of degree I [20 

r(q+Z+l) 
V. r(a+l) 



L ?(y) = tttStt $ H> « + 1; y) = ^ « + 1; y), [« = 0,1,2,..-] (131) 



which yield 

e 

T(n + e + B 2 ) 



Ui oc Ut = fl f^Z7^[ 2 ^ BL K* _1 (-^), ^ iV + ^ = ~* (132a) 



n=0 

and 



Ut oc f/r = r(n + " + 6 i) [2 ^ ]WLKB3 " 1(2 ^ ) ' if - ^ + f = (132b) 

n=0 ^ 
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The above types of finite series occur under the same conditions [that is, ±irj + (B 2 /2) = 
—£] as the ones obtained from the Baber-Hasse expansions in ascending power series of 
(z — Zq) [28j [3j [1]. However, the power series are more general, since the expansions in 
confluent hypergeometric functions require restrictions on the parameters of the CHE. On 
the other side, finite and infinite series solutions concerning the same set present different 
multiplicative factors, exp(iuz) or exp(— iuz), and so it is possible that only one of these 
solutions is bounded at z = 00. 

The expansions in Coulomb functions afford a second type of finite series. For instance, 

if B 2 ^ 0, -1, -2, • • • and B 2 + (Bi/zb) = -I in Ui, [£ = 0, 1, 2, ■••] (133) 

then the coefficient 7* of h\_ x is zero for n = £ + 1 and so the three solutions Ui are given 
by finite series with < n < £. For this second type of finite series, the hypergeometric 
functions do not reduce to polynomials and, in general, the solutions are not bounded for all 
values of z. Indeed, U\ is bounded for finite values of z, but may be unbounded at z = 00 
because the limit (1A.7|) for $ contains both the factors exp (iuz) and exp (—iuz). On the 



other side, the solutions Uf can be unbounded at z = because in general the functions \& 
present logarithmic terms at z = 0. 

The second type of finite series does not occur in solutions of the DCHE (zq — > 0) but is 
the only type which may be valid for the CHE with rj = and also in the Whittaker-Ince 
limit fl4]) of the CHE. This fact becomes relevant because these equations do not admit finite 
series in ascending powers of z or z — Zq. 

3.1.3 Solutions for 77 = and the spheroidal equation 

As in the case of two-sided solutions, if rj = we write the confluent hypergeometric 
functions of the first set in terms of Bessel functions, redefine the series coefficients as 

bl = i-T (n + f ) a\ (134) 

and add the expansions in series of Bessel functions of the second kind, Y K . Thus, from the 
first set f)122p we find the set given in Eq. (I137al) . The full group is generated by the 
transformations (TlBl applied as in Eqs. ( I114jl . that is, 

U ( i\ U { i ] = T^f, U { i } = T 2 U { i\ U { i ] = 7W ( J : first subgroup; 

= T 4 U { {\ U [ i ] = = T 2 U ( i\ U { {> = : sec. subgroup. 

Due to the redefinition f )134p . relations (I117P and (I118P give 
a ai + (3 a = 

a n a n+ i + (3 n a n + %a n _i = 0, [n > 1] J Pi- /3 2 - /3 3 - 



«07l «172 «273 / 10 _v 

Po = -o o o , (135) 



d ai + /3o«o = 0, 

• - 1 h ° 1 r° » 1 ° n L _v r «o [7i ~ a-i] «i72 
&ia 2 + Piai + [71 - a_ij a = 0, > p = « = , (13b) 



a n a n+ i + /3 n a n + %a n -i = 0, [n > 2] 



Pi- P~2~ 



where the coefficients a n and 7 n have also been redefined. For example, the redefinition 
(I134p leads to the redefinitions d* H- —i(n + B 2 /2)al and 7^ H- ij^/ (n — 1 + B 2 /2). This 
explains the difference between relations (11181) and (11361) when rj = and B 2 = 1. 
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Now we write first and the fifth sets of solutions. The first set reads 

oo 

^W=^£^a-iH. ( 137a ) 



n=0 



Mao[B+ i][n-%] . ,-H-B, .i = ^["+^-2]["+^-i+^] (137b) 

"n (2n+73 2 +l) ' ™ L' 4 ^ -°2 J- J -°3, | n (2ri+B 2 -3) ' 

For the fifth set, we find 

oo 

U%\z) = {z- z$-%^%?® x \uj{z - z )) , (138a) 



2 2 
n=0 



-^ [n+l][n+B 2 + ^] ° 5 . -a,s [n+B 2 -2] [n-1-^] /jggjx 

(2n+B 2 +l) ' ^ n (2n+B 2 -3) 

In both sets the independence of the Bessel functions requires that 

B 2 ^ 0, -1, -2, • • • if the series begin at n = in f/^and , (139a) 

while the recurrence relations are given by 

Eqs. (JUSD if 5 2 ^ 1, Eqs. ([136]) if B 2 = 1 for /7 ( ^and U {j 5 ] . (139b) 

On the other side, for the full group of solutions (Appendix C), we find that 

2ui + 2 ^ 0, -1, -2, ■ ■ - , [independence of in f/^and t/gil (140) 
for series beginning at n — 0, while the recurrence relations are given by 

Eqs. (II35D if 2^ + 2 ^ 1, Eqs. ([136]) if 2i/ t + 2 = 1, (141) 

the parameters being defined in (I116p . Relations (I139al) and (1139bj) stand for z/j = v\ in 
ffT40l) and ffTiil . 

For the present case (77 = 0), finite and infinite series are not possible in solutions be- 
longing to the same set, unlike the case 77 7^ 0. Further, finite-series solutions are of the type 
given in (11331) which may be unbounded at some of the singular points. This holds also for 
solutions of the spheroidal equation. 

In fact, in the literature we have found no finite-series solutions for the spheroidal equa- 
tion. The preceding solutions give four sets of solutions which do not admit finite series, and 
four sets which admit finite series, namely, 



only infinite series: S { {\ S%\ S%\ S^] 
finite or infinite series: S%\ S { {\ S%\ S { $ . 



(142) 



As examples we write the solutions S 1 and ■ Thus, by setting 00 = —27 in Eqs. (I58bj) 
the solutions (I58ap corresponding to read 



y-i 



^nlMv ~ !)]. \M ^ -1, "3/2, -2, • ■ ■ ] (143a) 



n=0 

where the recurrence relations for a\ are given by 

Eqs. (P35J if fi ^ -1/2, Eqs. (HH if fi = -1/2 (143b) 
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with 

^= 2J %XXT^ ft = (n + Ai)(n + M +l)-A, 7i= yg' ( 143c ) 

Due to the restrictions /i / -1, —3/2, —2, ■ • • which assure the independence of each term 
of the series, it is impossible to have 7* = and, consequently, these solutions do not admit 
finite series. Similarly, the other solutions which do not permit finite series are expansions 
in series of V'ndL anc ^ require restrictions on /i. 

The solutions admitting finite series are expansions in series of ipty which are all in- 
dependent, regardless the values of /i. In particular, from given in Appendix C, we 
find 

r 4- 1 1 £ 00 

(/*,*)= E a »^^-l)L (144a) 

-I n=0 

where the recurrence relations for are given by Eqs. f !135j) if the series begin at n = and 
by Eqs. (1128aj) if the series begin at n = N > 1. The coefficients for (Z^ 21X6 

^= 27( t':r' & = n(n+l)-A, 7* = ^Sf. (144b) 
So, if /x = — £ — 1 = — 1, —2. — 3, • • • , then 7^ +1 = implies finite series with < n < I. 

3.2 Solutions for the double-confluent Heun equation 

Only the solutions Ui, U3, U5 and U7 of the CHE admit the limit z ~~ * 0- Besides 
this, the limits of Ui and U3 are equivalent to the limits of U5 and U7. Thus, we obtain a 
subgroup of two sets of one-sided solutions for the DCHE. A second subgroup is generated 
by using the transformation t 2 . In fact, denoting by Ui the first set for the DCHE and using 
the transformations (I82a|) . the full group is given as in Eqs. ( 185|) . that is, 

UiCs), V 2 (z) = tiUi(z); \J 3 (z) = t 2 Ui(z), U 4 (z) = t 2 V 2 (z) = t 3 t 3 (z). (145) 

Conditions for connecting solutions of a given set may be found as in the case of the CHE. 
In the following we write the four sets of solutions and discuss the possibility of finite-series 
solutions (Heun polynomials). 

The limit z Q — > modifies the expressions for coefficients f ll23p but leaves the solutions 
(II 2 2 j) unaltered. So, we have 



00 

B 



t 1 (z) = z-^^2b 1 n %_ 1+Sz (r,,uz) & Eq.ro [B 2 ^ 0, -1. - 2. • • • ] (146) 



n=0 



or, explicitly, 



^T(2n + B 2 ) V 2 J 



n=0 

00 



(147) 



-r[nT^ + (B 2 /2)] 
where in the recurrence relations for b. 



v. 



□ 1 2«w_Bi(n+l) 

or ~ — 



'» (2n+B 2 )(2n+B 2 +l)' 



fa = n (n + £ 2 - 1) + B 3 + ( JXgg+k) > (148) 

ol _ 2iw_Bi[n+_B 2 -2][n-l+i»?+-^-][n-l-ir ? +^-] 
7n — (2n+B 2 -3)(2n+B 2 -2) 1 
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Thence, the recurrence relations are again given by (11261) . that is, 

7] ^ : Eqs. (HHD if B 2 ^ 1, 2; Eqs. ([USD if 5 2 = 1; Eqs. flTTSD if #2 = 2. 
77 = : Eqs. (dUD if 5 2 ^ 1; Eqs. dHHJ) if B 2 = 1. 

The second set of solutions, corresponding to the solutions U3 of the CHE, is 



(149) 



%{z) = tJJiiz) = J2 ° h l%+i-^ (V, uz) , [B 2 ^ 4, 5, 6, 



.E2. 

n=0 

which, explicitly, reads 



00 



U 2 = e*"*"**-* Yl r[L+rl> [" + 2 + ^ - f , 2n + 4 - B 2 ; -2zu;z] , 

n=0 

00 

^ = e^^* Yl r^T^] * + 2 ± ^ - f , 2n + 4 - 5 2 ; T 2iw*] 



(150) 



with 



o 2 _ 2i^_B! (n+1) 

— (2n+4-B 2 )(2n+5- J B 2 )' 



^ = -(n + l)(n + 2 - B 2 ) - 5 3 - ^Jl^Uy (151) 



^,2 



2iwBi [n+2-B 2 ] [n+l+i V - ^\ [n+1 -177- 



7n (2n+l-B 2 )(2n+2-B 2 ) 

in the following recurrence relations 

r] ^ : Eqs. ([117]) if B 2 ^ 2, 3; (UTS]) if B 2 = 3; dHHD if B 2 = 2; 
77 = : Eqs. (JTT7J) if £ 2 ^ 3, (HUD if fi 2 = 3. 



(152) 



If ?7 = 0, Ui and U 2 do not admit finite-series expansions (<j n 7^ 0) and can be rewritten 
in terms of Bessel functions. Finite series are possible if 77 7^ because, by setting v\ = 
(B 2 /2) — 1 and v 2 — 1 — (B 2 /2), we can interpret the solutions as in the two items following 
relations (11301) . 

For the second subgroup (U3 and U4) the transformation t 2 changes the form of the 
recurrence relations since ir] is transformed into (B 2 /2) — 1. Precisely, the relations (I117p . 
(I118P and (11191) are replaced by 

a k + f3 b = 0, a n b n+1 + j3 n b n + %b n ^ = 0, [n > 1] <S> flTT7J. (153a) 



1 + /3 &o = 0, 

aik + /§i6i + {71 + [(f " I)' - 1] ^0 = 0, (153b) 
a n 6 n+ i + /3 n 6 n + 7„6„_i = 0, [n>2]. 



a &i + 



a_i 



= 0, 

(153c) 

a„6 n+ i + $ n b n + 7„6„,_i = 0, [n > 1]. 
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If B 2 = 2, relations ( 1153d) reduce to ( I153ap . 
Up to a constant factor, the third set is 



e iuiz+^ z l- 



+ir) 



n=0 



B2 

2 



iB 1 
~2z~ 



where the conditions 2ir) 7^ 0, —1. — 2. 
begin at n — 0. Explicitly, 



[2i77^-2,-3,-4,...-] 
assure the independence of each term if the series 



U 3 (z) 



Ut(z) 



r[2n+U+2] ®l n + W 



n=0 



f , 2n + 2i V + 2; f ] , 



1 71 + 777+- 



(154) 



n=0 



* [n + z?7 + 1 ± 2n + 2ir] + 2; ±^l] , 

which correspond to solutions (186aj) with z/ = ^77. The recurrence relations for are 

B 2 ^ 2 : Eqs. (TI53aj) if 2zr/ ^-1,0, fTT53b]) if 2ir] = -1; (TT53c]l if 2ir/ = 0; 
5 2 = 2 : Eqs. ( TT53lj) if 2zr/ ^ -1; f fl53b|) if 2zr/ : 

with 

. 3 _ 2iw_Bi(n+l) 



1 



(155) 



a 



(2n+2iri+2)(2n+2ir]+3) ' 

/33 = 5 3+ (n + ^ + l-f) (n + iv + f) + {2 ^\Sl^) 



(156) 



7 3 



2iuBi [n+iri-1+^2-] [n+irj+l-^] [n+2irj\ 
(2n+2iri-l)(2n+2ir]) ' 



As in the case of two-sided solutions, the fourth set is obtained by replacing (77, uf) by 
(—77, —u) in (U 3 , U^) and respective recurrence relations, that is, 



U,(z),Ut(z) = t 3 U 3 (z),t 3 U*(z) 



[2^2,3,4. 



(157) 



These correspond to solutions (1HT1) with v = —it]. 

If B 2 = 2, U 3 and U4 do not admit finite-series expansions and can be expressed in terms 
of Bessel functions. Finite series are possible when B 2 7^ 2 since (for £ equals to zero or a 
natural number) 

• If irj+(B 2 /2) = —£, U 3 and U 3 are solutions linearly dependent and lead to finite series 
with < n < £, whereas the series in U 3 ~ begins at n — £ + 1. If 2 + irj — (B 2 /2) = —£, 
U 3 and U 3 ~ are multiple of each other and give finite series with < n < £, whereas 
the series in U 3 begins at n — £ + 1. 

• If —ir]+(B 2 /2) = —£, U4 and give finite series with < n < £, whereas the series in 
XJ~l begins at n — £+ 1. If 2 — ir] — (B 2 /2) = —£, U4 and are multiple of each other 
and give finite series with < n < £, whereas the series in begins at n = £ + 1. 

By Eq. ( 11311) . the finite series may be written in terms of generalized Laguerre polynomials. 

The convergence of the one-sided infinite-series solutions follows from the convergence of 
two-sided solutions (see section 2.C): 



Uf and Uf converge absolutely for \z\ > but may be unbounded at z = 00; U 3 and 
llf converge for \z \ < 00 but may be unbounded at z = 0. 
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• By the ratio test, the series which appear in Ui converge for any z. Nevertheless, it 
is necessary to compute the limits of U\ as z — > and z — > oo in order to get the 
behaviour of these solutions at the singular points z = and z = oo. 

Solutions for the Whitaker-Hill and Mathieu equations, regarded as DCHEs, may be 
obtained by using (|95|) and (11001) . On the other side, from the expressions for p^ and we 
see that U3 and U4 do not admit the Whittaker-Ince limit fl3]). Thus, in this limit we will 
find only the two sets written in Eqs. ( 1165al) and (I166al) . 



3.3 Solutions for the Whittaker-Ince limits of the CHE 

The Raabe test has been used in a previous paper [I] in the study of one-sided series 
solutions for the Whittaker-Ince limit of the CHE. Such solutions are exactly the ones ob- 
tained by applying the Whittaker-Ince limit ([3]) to the preceding one-sided solutions for the 
CHE. Thus, here we restrict ourselves to the form of the recurrence relations, write the first 
set of solutions and indicate how to obtain the other sets. We show as well how the domains 
of convergence of some solutions for the Mathieu equation are affected by the Raabe test. 

For series beginning at n — 0, instead of the forms (1117111191) for the recurrence relations 
we find, respectively, 

« ci + A)Co = 0, „ Sio'ji a 2 73 M ~o\ 

} Po = = g , (158) 

a n c n+ i + (3 n c n + 7„c n _i = 0, [n>l] J Pi- P 2 ~ Pz~ 



d ci + [3 c = 
&ic 2 + $\h + [ A 

+ %c n -\ = 0, [n> 2] 



- go [71 + a-i] «i7 2 a 2 7 3 
aic 2 + P1C1 + [71 + a-i] c = 0, } =>■ p = — , (159) 

PI— P2— P3~ 



a c± + 



$0 + a- 



c = 0, 



a 7i ana Qj 2 73 / 1Rm 
A) + «-i = -o o o , (160) 



+ /3 n C n + 7 n C n „x = 0, [n > 1] J /?3 

Once more, n > —1 for a n \ n > for /3 n ; and n > 1 for j n . If the series begin at n = N > 1 
(that is, if a N _i = 0), the recurrence relations have the form (I128aj) . 

Observe that the relations (I158H160]) are not obtained by applying directly the Whittaker- 
Ince limit on (111711119]) . For example, before taking the limits of the solution Ui given in 
(I122p . we redefine the series coefficients as - see Eqs. (I68ap and f!68bj) - 



Then, we insert these expressions into (111711119]) and perform the limits. 
The first set, corresponding to solutions (I122|) . is [1] 

oo 

U^(z) = z^~ B ^ ^(-^n^Bi-i (2v^) , & ^ 0, -1, -2, ■ ■ ■ ] (161a 

n=0 

where recurrence relations for c\, if the series begin at n — 0, are given by 



Eqs. (USE]) if B 2 ^ 1, 2; Eqs. if B 2 = 1; Eqs. (RM) if B 2 = 2, (161b) 
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with 



(n+l)L- 



(2n+B 2 )(2n+B 2 +l) ■ 



fa = n (n + B 2 

qz {n+B 2 -2)\n+B 2 



1) - 2|° 



qz (B 2 -2 



2{2n+B 2 -2){2n+B 2 )- 



From f> ( ^(z) 
First : 
Second : 



(2n+B 2 -3){2n+B 2 -2) ' 

the two subgroups are obtained by using the transformations fl66|) as 



(161c) 



u 



U) 



U) 



u 
u 



(j) 

2 



u 



U) 



3? 2 U%\ 



Mi) 



00 

6 ■ 



f7 



U) 
1 

(7) 



0'). 

3 ) 



(162) 



The convergence of the one-sided infinite series now results from the convergence (1741) for the 
two-sided series, the point z = oo requiring special attention. Explicitly: (i) the solutions 



U j in series of Bessel functions of first kind converge for finite values of z\ (ii) for j = 2, 3, 4 
the solutions U ^ converge for \z\ > \zq\ in the first subgroup and for \z — z \ > \zq\ in the 
second subgroup but, according to the Raabe test, they may converge also at \z\ = \zo\ and 

\Z - Zq\ 



\zq\ under conditions similar to ( 11291) . namely, 



\z\ > \zq\ if 



[Re 


B 2 H 




< 1 in 










1 Re 




_ Sx 


> 1 in 






z ° . 





rU) 



VCj). 

y 3,4> 



z \ > \z \ if 



fRe 




> 


— 1 in 




. z ° . 






1 Re 


§X 


< 


— 1 in 











V(i) 

y 5,8' 

V(i) 



where Um, n means: U m, and Un- The behavior of the solutions at z 



(163) 

oo must be determined 

by using the limits (1571) . In virtue of the multiplicative factors it is also necessary to examine 
the behavior of the solutions at the other singular points lying inside the regions afforded by 
the ratio test. 

Now we consider the convergence of the solutions w Mu) for the Mathieu equation. 



U;(z) where U;(z) are the solutions of the Whittaker-Ince 



According to ( 180]) . w \ \u) 
limit of the CHE with z Q = 1, B\ = —1/2, B 2 = 1 and z = cos 2 (au). Thence, the 
expansions of w^(u) in series Bessel functions of the first kind converge for all u but may be 
unbounded at the singular points, as state above. In comparison with the usual literature, 
the modifications implied by the Raabe test consist in the equal sign in the following domains 
of convergence (j = 2, 3, 4): 



cos(crti)| > 1 for w i and w { 



2 ■ 



sin <ju ) 



o (j) a 

> 1 for w 5 and w 



U) 



(164a) 



(i) °0') °(i) 



In effect, by taking a = i and writing explicitly the solutions, we find that w \' , w" 2 , w\ 
and w g are equivalent, respectively, to the solutions 28.23.6, 28.23.8, 28.23.7 and 28.23.11 
of Ref. [18]. For the other solutions, the conditions given in ( I163[) are not satisfied and, then, 



COS (JW 



> 1 for w 3 and w y A 



sin cm 



> 1 for w g and 



(164b) 



as usual. 



3.4 Solutions for the Whittaker-Ince limit of the DCHE 

The Whittaker-Ince limit of the one-sided solutions for the DCHE yields two sets of 
solutions for Eq. (jSJ), as mentioned at the end of section III.B. The first set can be obtained 
as limits of Eqs. (1161 alfTBTcj) when zq — > 0. Thus, 

oo 

U [ {\z) = z^/ 2 j^-ir^g^ (2V?S) , [B 2 + 0, -1, -2, • ■ •] (165a) 

n=0 
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where recurrence relations for c\ are given by 

Eqs. (USED if B 2 ^ 1, 2; ([159]) if B 2 = 1; ([1591) if B 2 = 2, (165b) 



with 



A, 1 — gffi(n+l) ol _ / , R _ \ , rj qB 1 (B 2 -2) 



(2n+B 2 )(2n+B 2 +l)' ^ n v 1 ^ ; 1 13 (2n+B 2 -2)(2n+B 2 ) ' 

^1 - gBi(n+B 2 -2) MfiVl 
'n — (2n+B 2 -3)(2n+B 2 -2)- l iUJ W 

The second set is obtained as U 2 {z) = TU^{z), where T is the transformation written in 
Eq. ffTUTl) . Hence 



?7 ( ?(z) = z^^/ 2 ^(-l) n c^^ +3 - B2 , [S a ^ 4, 5, 6, • • ■] (166a) 



r 0V^ - e Bl / z r(l-- B a)/2V^^_ 1 ^n«2 5 r(7) 

71=0 

where the coefficients c 2 satisfy the recurrence relations given in 

Eqs. (USED if B 2 ^ 2, 3; (H59J) if 5 2 = 3; jTBDj) if 5 2 = 2. (166b) 

with 

^,2 _ /i 2 - f« J. 1 Vn 4 9 - R U R gB 1 (B 2 -2) 

"n — (2n+4-B 2 )(2n+5-B 2 )' ~~ V 4 VV' t "r ^ T" -D3 (2n+2- J B 2 )(2n+4-B 2 ) ' 

= gBi(n+2-.B 2 ) / 16 g x 

In (2n+l-S 2 )(2n+2-B 2 )' l iUU W 

By examining the coefficients d n and j n , and taking into account the restrictions on B 2 , we 
see that the above solutions are given by infinite series beginning at n — 0. 

^From the convergence of the two-sided solutions for the DCHE, we find that 

• and U 2 converge for all z but, to decide if they are bounded or unbounded at 
the singular points, we have to compute the limits when z — > and z — > oo. 

• Ui and U 2 (j = 1, 2, 3) converge for \z\ > but their limits when z — > oo may be 
unbounded due to relations f[5Tl) . 



4 Possible Applications 

In this section we consider two examples which use solutions for the CHE. In the first 
example we discuss solutions for the Klein-Gordon equation for a scalar test field $ in the 
gravitational backgound of a singular and a non-singular spacetimes. As the time dependence 
of the scalar field obeys Mathieu equations without any arbitrary constant, we have to use 
two-side series solutions. The parameter v must be determined from the characteristic 
equation. Then, the Raabe test assures that the wavefunction for the non-sigular spacetime 
is bounded and convergent for all values of the time variable. 

The second example deals with the one-dimensional Schrodinger equation for a family 
of quasi-exactly solvable potentials. In addition to the expected solutions given by finite 
series, for a subfamily of the potentials we find infinite-series solutions which, due to the 
Raabe test, are bounded and convergent for all values of the independent variable. If the 
equation reduces to a Whittaker-Hill equation (WHE), we express such solutions by one- 
sided expansions in series of Coulomb functions; otherwise, the solutions are expressed by 
two-sided series. For the latter case the parameter v may be chosen conveniently, while the 
characteristic equation must be used for determining the energy levels. 
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4.1 Klein-Gordon equation in curved spacetimes 

In its conformally static form, the line element ds 2 = g liu dx^dx 1 ' for nonflat Friedmann- 
Robertson- Walker spacetimes is written as 



ds 2 = [A(r)f 



dr 2 — d\ l 



sin 2 (\/ex) 



(d9 2 + sin 2 6d<p 2 ) 



x»={t,xA<P) (167) 



where e = ±1 is the spatial curvature, r is the time variable, whereas x> $ an d V 9 are the 
spatial coordinates. The Klein-Gordon equation for a field $ with mass M (H — c = 1) is 

^ {V^gg^d,,®) + (m 2 + qR) $ = o, d„ = d/dx», 

where g is the determinant associated with g^, R is the Ricci scalar, g — 1/6 for conformal 
coupling, and £ = for minimal coupling. By performing the separation of variables 



*{x,0,<P,t) = [Air)]- 1 T(r) X( X ) 0(9) e 



imtp 



m = 0,±l,±2, ••• 



168) 



one finds that X and B are given by the same special functions for any scale factor A(r) 
j, while T obeys the equation 



dr 2 



+ 



k + M A + (6g - 1) 



1 d 2 A 



+ e 



T = 0, 



A dr 2 

The constant of separation k, determined from the spatial dependence of is given by 

k = 1, 2, 3, • • • if e = 1, and < k < oo if e = — 1. 

For a nonsingular model of universe and for (singular) radiation-dominated models, Eq. 
(I169P reduces to Mathieu equations. 

4.1.1 Nonsingular metric 

For the nonsingular case, the scale factor A{r) is given by [23] 

A(t) = a coshr, e = —1, — oo < r < oo, 
where a is a positive constant, leads to the modified Mathieu equation 



d 2 T 
~d^ 

So, in Eq. (j6]) we have 
a = i, 



\m 2 &\ + ^M 2 a^ cosh(2r) 

_ _ 



T = 0. 



-k 2 - (M 2 a^/2) , k = Ma /2, 



U = T. 



Solutions for this problem have already been given in Ref. [25] where the convergence at 
r = is not discussed. Here this question is solved by using the Raabe test. 

From the Heine-type solutions, only iy y given in equation (18 lap afford convergent and 
bounded wave functions for all r G (— oo, oo). We find the solutions 



T {j) (r) = (-iAn^lU+i (Maocoshr) , [2u £ Z] 



(169a) 
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where the recurrence relations for c n are 



M 2 &1 c n+1 + [(An + Au + 2) 2 + Ak 2 + 2M 2 a 2 ] c n + M 2 a 2 c n _i = 0. 



(169b) 



The relations among the Bessel functions [H] imply that only two of the four solutions (I169ajl 
are linearly independent. Similar results are found by treating the Mathieu equation as a 
CHE. In effect, by using given in (164aj) we obtain 



T0(t) = a nZ^l +{1/2) [2Ma cosh 2 (r/2)] , [2u £ Z] 



(170a) 



where the recurrence relations for 

M& [n + v + l]a n+ i - (n + i/ + |) 2 + k 2 + M 2 a\ a n + Ma [n + z/]a n _i = 0. (170b) 

Thus, the Raabe test assures that and T^, as well as the corresponding wavefunctions 
f)168p . are bounded and convergent for all r G (—00,00). Notice that the solutions fllOlal) . 
resulting from the DCHE, do not satisfy the conditions fll01c|) for all r G (—00, 00). 

4.1.2 Singular metric 

For radiation-dominated spacetimes, A{r) = a sin (er)/ y^e (r > 0) and, so, 

j| + [« 2 + |M 2 a 2 - ^M 2 a 2 cos (2^r)] T = 0. (171) 

We consider only the case e = —1. We take a = z, a = (M 2 ap/2) — k 2 , k = Mao/2 and 
u = r. Then the Heine solutions u; ^ (it) given in equation ( 18 lap lead 



r 0-) (r) = ^ (-l) n c n Z^ +2 , + i (Maocoshr) , \2v £ Z] (172a) 

n=— 00 

where the recurrence relations are 

M 2 a^ c n+ i+ [(4n + Av + 2) 2 + 4k 2 - 2M 2 a^] c n + M 2 a\ c n _i = 0. (172b) 

On the other side, from the solutions w 1 (it) given in ( I64al) (CHE) we find 

00 

T®(t) = Yl a «^i+(i/2) [2Ma cosh 2 (r/2)] , [2v £ Z] (173a) 

n=— 00 

where the relations for 

Ma [n + v + l]a n+ i - (n + z/ + ±) 2 + K 2 a n + Ma [rt + i/]a n _i = 0. (173b) 

Once more, T^'> and T^ are convergent and bounded for all r > but now the wavefunctions 
(11681) become unbounded at r = due to the factor 1/A(r) = l/(aosinhr). Again, the 
solutions (1101 al) resulting from the DCHE are inappropriate for e = — 1. 

Therefore, if e = — 1 the solutions for the modified Mathieu obtained from the CHE and its 
Whittaker-Hill limit are suitable for the singular and nonsingular metrics, while the solutions 
obtained from the DCHE are unsuitable. For the singular metric, the unboundedness of the 
solutions fl!68p at r = is expected since at this point there is a physical singularity in the 
sense that the pressure and density energy diverge. 
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4.2 Schrodinger equation for quasi-exactly solvable potentials 

Now we consider problems involving solutions given by finite and infinite series for the 
CHE. For this end, we write the one-dimensional stationary Schrodinger equation for a 
particle with mass M and energy E as 

g + [£ - V(«)ty = 0, u = ax, £ = |^£, V(u) = |^V(x), (174) 

where a is a constant with inverse-length dimension, h is the Plank constant divided by 2ir, 
x is the spatial coordinate and V(x) is the potential. For V(u) we choose the Ushveridze 
quasi-exact solvable potential p2] 



V(u) = 4/3 2 sinh 4 u + 4/3 [/3 - 2( 7 + 5) - 2£] sinh 2 u + 4 [6 - J] [6 - §] 



sinh « 



4 [7-J][7-|]^, [i = 0,l,2,...] (175) 

where /3, 7 and 5 are real constants with j3 > and 5 > 1/4. 

When 5 > 1/4 and i is zero or a natural number, the above family of potentials is 
quasi-exactly solvable because it admits bounded wavefunctions given by finite series which 
allow to determine only a finite number of energy levels. However, for 1/4 < 5 < 1/2 and 
l/2<5<3/4we also find infinite-series solutions which are convergent and bounded for all 
values of the independent variable: this suggests the possibility of determining the remaining 
part of the energy spectra as solutions of a characteristic equation. For 5 > 3/4 we find no 
solutions like these. 

Notice that Ushveridze supposed that I = 0, 1, 2, ■ • • , but we will find that if 

(7, = (U), (L»> (I.I). (176) 

then we also have quasisolvability when I is a positive half-integer. In addition, Ushveridze 
supposed that u £ (—00, 00). However, we get 



lim V(u) = 00, lim'V(u) 

u— >±oo u— >0 



3. 

•i' 



-4 [7- J] [7-!], if5 = ior5 
-00, if*e (177) 
[+00, if ^ [If]. 



Hence, for 5 ^ [1/4,3/4] there is an infinite barrier at u = and, so, we can suppose that 
u > or u < 0. 

4.2.1 Wavefunctions for the Whittaker-Hill equation 

If 7 and S take the values (11761) . the potential reads 



V(u) = 4/3 2 sinh 4 u + 4/3 [/3 - 2(7 + 5) - 2£]sinh 



2 «, 



°' 2' X ' 2'- 



(178a) 



where u £ (—00, 00). Thence, by using sinh 2 w = [cosh(2w) — l]/2 and sinh 4 w = [cosh(4w) — 
4cosh(2w) + 3]/8, Eq. fl 174ft becomes a modified WHE ((Jj) with the parameters 

a = i, ■&= _ £ + 4/ 3(£ + 7 + 5 ) ; p+ i = 2 (£ + 7 + 5), £ = 2/3. (178b) 

In fact, the WHE also occurs in the cases of the Razavy potential [26J and the symmetric 
double-Morse potential considered by Zaslavskii and Ulyanov 
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On the other side, the substitutions 

^ = cosh 2 w, ip(u) = tp{u(z)} = U(z), [z>l] (179) 

transform the Schrodinger equation for the preceding potential into the CHE (CQ) with 

zo = l, #i = -§, B 2 = l, £? 3 = f, kj = ±0, ir 1 = ±{t + 1 + 5), (180) 

where the plus or minus sign must be used throughout. Thus, we can attempt to solve 
the problem by using known solutions for the CHE. For example, from the Baber-Hasse 
expansions in power series, the Hylleraas solutions or from the Jaffe solutions [U 1221 ED] 
we obtain even and odd finite-series solutions bounded for z > 1: such solutions allow to find 
only a finite number of energy levels. There are also infinite-series solutions which, however, 
must be discarded because they are not bounded for any admissible value of z. 

On the other side, if we use one-sided series solutions tSi(z) in terms of Coulomb wave- 
functions, we shall find that (for £ equal to a non-negative integer or a positive half-integer) 

• Even and odd finite-series solutions which are convergent and bounded for all z > 1; 
these finite-series can be written in terms of generalized Laguerre polynomials by means 

of (031). 

• Even infinite-series solutions which, due to the Raabe test, are convergent and bounded 
for all z > 1. 

• Odd infinite-series solutions which converge and are bounded only for z > 1, and odd 
solutions which converge and are bounded only for finite values of z. To cover the 
entire interval z > 1 it is necessary to consider two of such solutions. 

Therefore, in this case we could find additional energy levels by solving a transcendental 
characteristic equation. 

In effect, since B\ = —1/2 and B 2 = 1, we find that all the confluent hypergeometric 
functions which appear in Uj(z) are linearly independent. It is sufficient to consider the 
solutions Uj = (Ui, U^ 1 ) with % — 1, 2, 3, 4. We take 7 = 5 = 1/4 and consider the other cases 
in the paragraph containing Eqs. 0186]) . As the expansions IJi(z) are not affected when we 
change (07,77) by (— o>, —77), in Eqs. (I180p we select 

iw = -0, if] = -I- (1/2), [ 7 = 5=l/4]. (181) 

Then, by using Ui = (JJ\, f/j 4 ) given in (I122p . we find 

00 

Mu) = e-^ 2 ^ ^g^ $[n-£,2n+l;2/?co S h^] , 

n=0 (182) 
#(«) = Y, rtgg^j] * [n + l^(£ + D ,2n+ 1; ±2/3 cosh 2 „] 

n=0 

which are even solutions with respect to u — > —u. According to Eqs. (11231) and f)126p . if the 
series begin at n = the coefficients b\ satisfy the recurrence relations (11181) with 

= K = -n 2 -\ +P {£+$), Y n = l(n + £)(n-£-l) 

If i is a non-negative integer ipi cx ipf represents bounded finite-series solutions with < 
n < £, while tp± represents infinite-series solutions with n > £ + 1 which become unbounded 
when z = cosh 2 w — > 00 due to the factor exp (f3z) {(5 > 0). On the other hand, if £ is 
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a positive half-integer represents infinite-series solutions convergent for z > 1 because 
Re[i?2 + (Bi/zq)] = 1/2 < 1 as required by the Raabe test; these solutions are bounded also 
at z = oo since relation (1A.6|) leads to 

bl „ .13 



n=0 L J 

By using ^(z) given in (IC.lj) . we again get even solutions but now the finite series hold if I 
is half-integer whereas the infinite series hold if £ is integer. We find 

oo 

^2 = e-P ^ « cosh it £ [n - £ + |, 2n + 2; 2/3 cosh 2 u] , 

n= ° (183) 

4± = e ^ cosh2u coshw V M 2Pc T % tf [n + 1 T U + I) , 2n + 2; ±2/3 cosh 2 «] , 

ra=0 

where, for series beginning at n — 0, the coefficients fe 2 satisfy the relations (I119p with 

2 _ _P f>2 _ _ ^ _l II 2 _ £ _l R \0 ^ II 



f> = — 



Thus, if £ is a half-integer then ^2 oc gives bounded finite-series solutions with < 
n < £ — (1/2). If £ is integer ^1 gives infinite-series solutions convergent for z > 1 because 
Re[i? 2 + (-Bi/^o)] < 1 as required by the Raabe test; these solutions are bounded also at 
z = 00 due to relation (1A.6|) . 

Odd solutions are generated from the expansions U3 and U4, but in these cases we find no 
single infinite-series solution converging for all z > 1 because the Raabe test is inapplicable 
for Uf- and Uf . Thus, from U 3 given in (lC.3p . we find the odd solutions 

00 

^ = e -0 cosh^ u sinh {2u) J- [ n + 1 _ ti 2n + 3; 2/3 cosh 2 u] , 

n=0 

*^-'«wi; ^* (184) 



^[ n + | T (£+i) ; 2n + 3; ±2/3 cosh 2 u] , 

where, for series beginning at n — 0, the 6^ satisfy the relations f 1 1 1 7 j) with 

= "/> 10by ^ = -(n ± l) 2 - § ± /3 (£ ± 1) , 7^ = . 

So, if £ is integer then ^3 oc ^3 is given finite-series with < n < £ — 1 (£ > 1), whereas ifjj 
is given by infinite series with n > £ (unbounded at z — 00). If £ is half-integer the three 
solutions are given by infinite series and ^3 stands for z 7^ 00, i/j^ stands for z 7^ 1 whereas 
ipz{u) (linear combination of "03 and i/j^) stands for z 7^ 1 and z 7^ 00. 

By using tl^z) given in C IO. 5 j) . we again get odd solutions but now the finite series hold 
if £ is half-integer whereas the infinite series hold if £ is integer. We find 

00 

^ = e -^ osh2 « sinh u hi [ "r[2n+ 2 ] 2n] " <l> - £ ± |, 2n ± 2; 2/3 cosh 2 u] , 

n=0 (185) 
ft = eT , cosh* « sinh u J- |gfgf^| * [n + 1 =F (* + |) , 2n + 2; ±2/3 cosh 2 u] , 

n=0 
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where the coefficients b\ satisfy the recurrence relations (11191) with 

- 4 n = -I im\ > & = & 7* = f [n + € + J] [n - £ - |] . 

Thus, if £ is integer it is necessary to take two infinite-series solutions to cover the interval 
z = cosh 2 w > 1 since ip±(u) holds for z ^ oo, ip~l(u) for z ^ 1, while ^j(w) holds for z ^1 
and z ^ oo. 

For the other values (11761) for the pair (5,7), the solutions can be obtained from the 
above ones (7 = 5 = 1/4) by a change in the parameter £. In effect, we find that 

(7, 5) = (i, f) or (|, 1) = -I - 1; (7, 5) = (f , f ) mj = - |, (186) 

without any other change in the parameters (I180p . Hence, by comparing the previous ex- 
pressions for it] with the one given in (I18ip . we see that in the above solutions it is sufficient 
to replace: (i) £ by £ + (1/2) in the first two cases, (ii) £ by £ + 1 in the last case. Anyway, the 
important conclusion is the same: there are even infinite-series solutions that are convergent 
and bounded for all z > 1. 



4.2.2 Wavefunctions for the cases 1/4 < 5 < 1/2 and 1/2 < 5 < 3/4 

For the Ushveridze potential (I175p . the substitutions 

^ = cosh 2 w, ip(u) = ip[u(z)} = z"!-*{z - l) s ~ l iU(z), [z>l] (187) 
transform the Schrodinger equation (j!74p into a confluent Heun equation with 



Zq = 1, Bi = -2 7 , B 2 = 2 1 + 25, B 3 = § + (7 + 5 - ±) 



2 



zw = ±/3, ir] = ±{£ + 5 + ~f). (188) 
Now we exclude the cases (11761) and suppose that £ is a non negative integer. We select 

iu = —(3, irj = —£ — 7 — 5. 

Then, by using for U(z) the power series written in Eqs. (29a-b) of Ref. [3], we find 

1 

^ hcT [u(z)} = e-P'zi-^z - l) s -*J2*n(z - l) n , [^ = 0, 1,2,-.-] (189a) 



n=0 



a n + 



where the series coefficients satisfy (a_i = 0) 

(n + l)(n + 25)a n+1 + \n(n + 2 7 + 25 - 1 - 2/3) + f + (7 + 5 - \) 2 - 2/35 

-2/3(n-£-l)a»_i = 0. (189b) 

Since /3 > 0, the previous finite-series solutions are bounbed for all z > 1 provided that 
<5 > 1/4. According to theory of three-term recurrence relations [7j, ^J 3aber represents £ + 1 
distinct solutions, each one with a different energy. 

By the other side, we find cases for which there are infinite-series solutions appropriate 
for any z > 1. For this we insert into (I187p the two-sided solutions and U 2 + given in 
(jl~9j) and (j23|) . respectively, and use the Raabe test along with the limit (1A.6I) . Thence by 
scribing convenient values to the parameter u, we find the solutions ipf and ip 2 having the 
following properties 
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• The solutions ipf are convergent and bounded for all z > 1 ifl/4<5<l/2. 

• The solutions ip2 are convergent and bounded for all z > 1 if 1/2 < 8 < 3/4. 

• If 5 = 1/2, then ifjf = ip^. For this case the Raabe test is inconclusive as to the 
convergence at z — 1. 

In effect, if U(z) = U^(z) we obtain 

K Wz] n 



oo 



^+( z ) = e -^z u - s+ ^(z-l) s -^ ^ 



n= — oo 



T[n + v + 1 + i + 7 + 5} 
x V[n + v+l -£-j-8,2n + 2u + 2;2^z], 1/4 < 5 < 1/2 (190a) 

where, in the recurrence relations f JT5|) for b\ we have 

1 2/3[n+t/+2-7-<5][n+i/+l+7-5] 
a n (2n+2i/+2)(2rc+2i/+3) ' 

= 'I + W + 7 + 6) - [n + „ + §] 2 - ^^K;^ ' ( 190b ) 

1 _ 2/3[n+i^-l+7+(5][n+^-7+(5][n+^+£+7+<S][n+i/-£-7-5] 
~ [2n+2u-l][2n+2u] ' 

By the Raabe test the condition 5 < 1/2 assures that the series converge at z — 1, while the 
condition 5 > 1/4 assures that the factor (z — l) 5 "^ 1 / 4 ) is bounded at z — 1. Similarly, if 
£/(z) = U£{z) we obtain 



r[n + z/ + l+£ + 7 + 5] 
x *[n + i/+l-^-7-5,2n + 2i/ + 2;2/3z], l/2<<5<4/3 (191a) 

where, in the recurrence relations ( 1151) for 65? 



2 _ 2/3[n+i/+7+3] [n+i/+l- 7 +6] o2 _ al 

"n (2n+2i/+2)(2n+2!/+3) ' , 

2 _ 2fi[n+^+l-7-<5] [rt+t/+7-<5] [ra+t/+l+7+<5] [n+t/-l~7-<5] ^ ' 

%i ~ [2n+2i/-l] [2ti+2i/] ' 

where = is a functional identity; in fact, and are different of each other because 
they hold for distinct intevals of 5. 

To assure that all the terms of the series are linearly independent and that the summation 
extends from minus to plus infinity, the parameter v must be chosen such that 

2u, v ± (7 + 8) and v ± (7 — 8) are not integers, (192) 

where the values for 5 are different for solutions (I190ap and (1191 al) . The linear independence 
is assured by requiring that 2v is not integer, without any restrictions on the parameters 
of the potential. Thus, for fixed values of 7 and 5, we can choose for v any value in the 
open interval (0, 1/2) convenient to satisfy the above conditions. The use of one-sided series 
would lead to restrictions on 7 and 5. 



5 Conclusion 

We have dealt with the convergence of Leaver's expansions in series of Coulomb wave 
functions for solutions of the CHE and DCHE. By redefining the Coulomb functions, we 
have completed the proof of convergence delineated by Leaver and, in addition, have found 
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that the Raabe test improves the regions of convergence for solutions of the CHE and its 
Whittaker-Ince limit (j4]) if certain conditions are fulfilled. It is worth noticing that in using 
the convergence tests we suppose that the independent variable z is finite. So, when z tends 
to infinite the behavior of each solution must be analyzed carefully. 

Other points that have also extended the meaning of Leaver solutions are: (i) the trans- 
formations of variables which generate new solutions (Leaver has used only transformation 
t2 for the DCHE); (ii) the Whittaker-Ince limit which gives solutions for Eqs. fll]) and ([5]); 
(iii) the construction of one-sided series solutions which, in most cases, lead to finite series 
solutions. These issues have been partially discussed elsewhere, but now they incorporate 
the consequences of the Raabe test for the solutions of the CHEs. 

The transformations of variables lead to solutions with different domains of convergence 
and/or different behaviors at the singular points. By this procedure, from the solutions 
for the Whittaker-Hill limit of the CHE we can recover all the Heine solutions in series of 
Bessel functions for the Mathieu equation [18] but the Raabe test improves the regions of 
convergence for some of these solutions. For example, in section 4, we have found that 
the Klein-Gordon equation in a non-singular model of universe has solutions bounded and 
convergent for all values of the time variable. 

In section 4 we have also considered the Schrodinger equation for the QES potentials given 
in Eq. (I175p . For certain intervals of the parameter S we have found some eigenfunctions 
given by one-sided and two-sided infinite series which, by the Raabe test, are convergent and 
bounded for all values of the independent variable. This type of solutions are not obtained by 
using other known solutions for the CHE as, for example, the Hylleraas and Jaffe solutions 

m Ea eu]. 

In fact, there is another QES problem [I] admitting well-behaved infinite-series eigen- 
functions, and these result from solutions for the Whittaker-Ince limit of the CHE. However, 
for the problem of section 4, the expansions in series of Coulomb functions do not give ap- 
propriate infinite-series wavefunctions if S > 3/4. We do not know whether it is possible to 
find new solutions for the CHE suitable for this case. 

On the other side, as particular cases of the solutions for the CHE we have obtained 
the Meixner two-sided solutions in series of Bessel functions for the spheroidal equation [9]. 
These solutions give one-sided infinite-series solutions which, in turn, lead to finite-series 
solutions as discussed after Eqs. (I142p . As far as we are aware, up to now these are the only 
solutions for the spheroidal equation which admit finite series despite the fact that we have 
found no application for them. 

We have also seen that the Raabe test is useless for solutions of the DCHE. On the other 
hand, we have found two subgroups of solutions for the DCHE - see Eqs. ( 18"5]) and ( I145p . 
One subgroup is obtained from solutions of the CHE when z — > 0; the other follows from 
this subgroup by the transformation t 2 given in Eqs. (I82ap but cannot be derived as limit 
of expansions in series of Coulomb functions. Thus, we can ask if there are solutions for the 
CHE which yield such subgroup for the DCHE when zq — > 0. 

It seems that is possible to solve the preceding problem by further investigation on known 
solutions in series of Gauss hypergeometric functions for the CHE [TU[ [3U [32]. For instance, 
we have the two-sided series expansion (TU] 

U+(z) = e^J2 < F {f -n-u-l,n + v + f ; #2 + f ;1- , (193) 

n 

where F(a, b; c; () denotes the Gauss hypergeometic functions, and the coefficients d\ are 
multiples of the coefficients b\ which appear in solutions ffT9l) . This solution is convergent 
for finite values of z and, so, cover a region of the complex plane left off by the expansions 
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in series of Coulomb wave functions. On the other side, from the limit [20] 

lim F [a, b-c-l- (c/y)} = y b ^(b, l + b-a;y), 

c— >oo 

we find lim 20 ^ Uf(z) oc U^(z), where U^(z) is just the solution given in Eqs. (I86ap for the 
DCHE, which does not result of any expansion in series of Coulomb functions. To get the 
solutions U?,(z) and U^{z) for the DCHE, we would need to find two additional expansions 
in series of hypergeometric functions for the CHE. We could also examine if such solutions 
supply the Meixner expansions in series of Legendre functions for the spheroidal equation 
[9J , having the same coefficients as his expansions in series of Bessel functions. 

Finally, notice that we have not considered the second type of Leaver solutions in series of 
confluent hypergeometric functions for the CHE. In fact, these solutions have been discussed 
in a previous paper [3] where, however, the Raabe test has not taken into account. First, 
we can show that the test is applicable to the expansions in series of irregular confluent 
hypergeometric functions. Second, for the DCHE the transformation £ 2 gives solutions which 
cannot been derived from Leaver's solutions of the CHE through the limit z — > 0. Thence, 
we have another unsolved problem connected with Leaver's solutions. 



A Confluent-hypergeometric and Coulomb Functions 

Here we write some useful formulas concerning the confluent hypergeometric functions 
and, in Eqs. fl A. 9j) and flA.101) . redefine the Coulomb wave functions. At the end we obtain 
the relations (1A.18|) which are important to apply the convergence tests for infinite-series 
solutions of the CHE and DCHE. 

The regular and irregular confluent hypergeometric functions, <3>(a, c;u) and \l/(a, c;u), 
are solutions of the confluent hypergeometric equation 



V& + {C-V)%-°<P = 0. (A.l) 

The functions <3>(a, c;y) and ty(a,c;y) are also denoted by M(a,c,y) and U(a,c,y), respec- 
tively [T7j. In fact, the following four solutions for Eq. ( 1A.1I) 



<pW(y) = $(a, c; y), <pW(y) = (a, c; y), 

(A.2) 

tpWty) = e »W(c - a, c; -y), ^\y) = y l ^{\ + a - c, 2 - c; y), 

are all of them defined and distinct only if c is not an integer [2D]. Different forms for </?W 
follow from the Kummer transformations 

$>(a,c;y) = e y $(c-a,c;-y), *(a, c; y) = ^"^(l + a - c, 2 - c; y). (A.3) 

In this article, we use only tp^, y^ 2 * 1 and p( 3 \ Their Wronskians are [20] 



W [<pV>, <p&] = W Ma, c; y), 9(a, c; y)] = -fg y~ c e\ 

W [</? (1) , <^ 3 )] = W [$(a, c; y), e y ^(c - a, c; e ±in y)} = r g^ y e Tinc y' c e y , (A.4) 

W [</? (2) , tpW] = W [*(a, c; y), e y ^(c - a, c; e ±in y)} = e ±in ^ y~ c e y . 

Therefore, if a, c and c—a are not zero or negative integers, the three solutions are well defined 
and any two of them form a fundamental system of solutions for confluent hypergeometric 
equation [20]. If c is not zero or a negative integer, the solutions are connected by [2D] . 

^ = ^ 9(a, c; y) + 9 (c - a, c; e^y) , (A.5) 
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When y — > oo, the behavior of of \I/(a, c; y) is given by [IT] 



y(a,c;y)~y > 

^— ' m! 

m=0 



37T 37T 

y < argy < — ; 



(A.6) 



while the behaviour of $(a, c; y) is given by 



$(a,c;j/) e y y a c ^ (1 - a) m (c- o) m _ 



T(c-a) ^m=0 



E 



r(c) 

oo (a)m(a— c+l) m 



r a 



m=0 



m! 



-y)" m , a^O, , c-a^0,-l, 



(A.7) 



where the upper sign holds for — 7r/2 < argy < 3ir/2 and the lower sign, for — 3ir/2 < argy < 
— 7r/2. In these limits (x) m denotes the Pochhammer symbol whose definition is 



(x) = l, (x) 1= x, (x) m = x(x-l)(x-2)---(x + m-l) =r(x + m)/r(x). 
On the other hand the Coulomb wave functions are solutions of the equation 



dy 2 + 



1 _ £1 

y 



2ri {n+v){n+v+l) 



0. 



(A.{ 



If r\ = 0, this can be written in the usual form of the Bessel equation by a substitution of 
variable. If rj ^ the solutions & n+v (y) = & n+u (r),y) are written in terms of one regular 
confluent hypergeo metric function $ and two irregular functions that is, 



where, by definition, we take 

<f>n+»(v, V) = r[ 2n+L+2] [^y] n+u+1( ^[n + v + l + ir 1 ,2n + 2v + 2- -2iy], 
±2ie r?7r e ±iy 

^n+v^^y) 



(A.9) 



T [n + v + 1 =)= irj\ 



[-2iy] n+u+1 ^/[n + v + l±ir],2n + 2v + 2; ^2iy\. 



(A.10) 



In V'n+i/ the irrelevant factors ±2zexp(?77r) are maintained just to connect the above defini- 
tions with the ones used by Leaver. In fact, for % n+u Leaver used the functions F n+U (r),y) 
and G n+U (rj,y), defined as 



F n +u(v,y) 



[T(n+u+l+iri)Y{n+v+l-i-q)] 1 / 2 
2e 7r i/ 2 r(2n+2^+2) 



e iy (2y) n+u+1 $(n + v + 1 + ir], 2n + 2v + 2; -2iy), 



(A.n; 



Thus, 



G n+v ( V ,y) ±iF n+u ( V ,y) = e^ 2 e 



tttj/2 =F«i"(n+f+l/2) 



r(n+t/+l±»?7) 
r(n+i/+l=Fw;) 



,1/2 



X 



e ±iy (2y) n+v+1 m [n + v + 1 ± 177, 2n + 2z/ + 2; ^2?y] . 



(A.12) 



and are obtained by dividing the above expressions by T n , defined as 



T„ = (l/2)e~ rin/2 (-i) n+u+1 [T(n + v + 1 + irj)T(n + u + l- ir])} l/2 . (A.13) 
Inversely, when n+ „ and ip^+v are multiplied by T n , we recover the Leaver normalization. 
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^From the properties of the functions F n+U (r],y) and G n+U (r],y) given in Eqs. (126) and 
(125) of Leaver's paper, we find that the functions (1A.9j) satisfy the equations 



d% n+v i(n + y){n + v + 1 + ir))(n + v + 1 — ir]) 



dy 



(n + ^ + l)(2n + 2z/+l) 



+ 1 



(n + f)(n + v + 1) 



i(n + v + 1) 



[n + u)(2n + 2u + 1] 



(A.14) 



and 



(n+u)(n+v+l+iri)(n+u+l— irj) <q, 
(2n+2u+l) 



{n+v)(n+v+l) 

y 



V 



(2n+2v+l) 



^n+v-l = 0. 



(A.15) 



For v and r\ fixed, by dividing all terms of (1A.15I) by (n 2 /2)% n+u and letting n — > ±oo we 
find 



whose solutions are 
and 



2i [1 _ _L1 ^ ^ j/. \i , JL1 

j/ L 2n\ w 2i L ~ 2n\ 



(A.16) 
(A.17) 



provided that y/n 2 = when n — > ±oo (this condition is satisfied if y is finite). Thus, there 
are two possibilities for the ratios between successive Coulomb functions. By demanding 
that these relations are valid also for 7] = 0, we find only one ratio: (i) the first expressions 
in flA.16j) and (I A . 1 7[) hold, respectively, for 4> n+v and ipn+v when n — > oo, (ii) the second 
expression in (1A.16|) is valid for the three functions when n — > — oo. In other words, 



q>n+v 



2n 2 



[i-H 2 "+§)] 



n + v + 1 2i 



2in z 



[1 + I (2v + I)] , W n+V = (0 n+ „ ^t +v ) , [n -> -oo]. 



(A.18) 



The above conclusions are obtained as follows. In the first place, if r] = the functions 
') n+u and V'n+i/ can be rewritten in terms of Bessel functions since 



$(n + i/ + 1, 2n + 2v + 2; -2zj/) = T [n + z/ + (3/2)] fe/^p-"-* e~^J n+u+ i (y), 
9(n + v + l,2n + 2v + 2; -2iy) = 



^e- i y + ™<^**\2y)-™-lHW v+1 (y), (A.19) 



tf(ra + v + 1, 2n + 2v + 2; +2i?/) 



V{2yy^-iHflJy) 



where J K is the Bessel function of the first kind, and Hjp and are the first and the 
second Hankel functions. Thence 

</w(o,y) = T^VyJn+^iy), i&„(o,y) = r^VvKit^ ( A - 20 ) 

where the constants C and C ± do not depend on n. In the second place, if y is bounded and 
k — v oo [7] 



(2), 



(A.21) 
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Combining ( 1A.20I) with ( 1A.21I) . we establish ( 1A.18I) for k = n + v + (1/2) when n — > oo 
(?7 = 0). On the other side, if ac — >■ — oo, we use the previous relations for H K ' (y) in 
conjunction with [18] 



HV(y) = e**"H}?'(v), H™{y) 



r(2), 



(A.22) 



Thus, we find ( 1A.18I) for ^j+j, = t/Wj, when k = n + v + (1/2) with rt — >■ — oo (77 = 0). For 
^ n+ ^ = 4> n+u , if y is bounded and k — > —00, once more we use the relation given in (1A.21j) 
for J K (y) since 



My) 



m=0 



m ! T(K + m + 1) \2/ 



2/\ 



2?fl 



r(« + i) 



m=l 



m ! T(K + m + 1) V2/ 



2A 



2 m 



(A.23) 



In this manner, we establish the ratio (1A.18|) for the three Coulomb functions when n — >■ —00. 



B Recurrence Relations for the Series Coefficients 



Now we present the derivation of the recurrence relations for the series coefficients of 
the two-sided solutions Ui(z) of the CHE, as well as of the one-sided solutions Ui(z). The 
relations for the other sets of solutions may be obtained from these by transformations of 
variables. Notice that: (i) the derivation is formal in the sense that, in each series, we 
suppose linear independence of all Coulomb wave functions; (ii) the importance of the three 
types of relations for one-sided solutions becomes apparent in the cases of the Whittaker-Hill 
and Mathieu equations. Details are afforded throughout the paper. 

The Leaver substitutions [1] 



U(z) = z~ B * /2 H(y) 

transform the CHE flTJ into 

y(y - ouz ) [0 + ( 1 

Ci = B\ + B 2 z , 
Expanding H(y) as 

00 

H(y)= b l n W n+u fay) & V x {z 

n=— oo 

and using Eqs. flA~8|) . (TOij) and fTA~T5|) we find 



y 



uz 



C 2 + ^ 

1 y 



H = where 



C 2 = B 3 - f [f - 1] , G 



B2Z 



Mi 
2 



Mi 



Z 2 



fay) 



Yl b n [ a l-l^n+u-lfa 2/) + Pl%+vfa y) + Ifl+lK+.u+lfa 2/)] = 0, 



(B.l) 



(B.2) 



(B-3) 



(B-4) 



where otn \ Pn^ and 7« are defined in Eqs. fl20|) . 

If v is such that the summation runs from minus to plus infinity, the preceding equation 
takes the form 



J2 [a l n b l n+l + (3lb l n + T * bl_ x ] W n+u ( V , y) = 0, 



(B-5) 



n=— oo 
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which is satisfied by the three-term recurrence relations fll5p provided that all the functions 
^ n+i ,(?7, y) are linearly independent. 

On the other hand, the series summation can begin at n = if v is chosen conveniently. 
In this case Eq. ( 1B.4|) can be rewritten as 



m=— 1 



m=0 



m=l 



or, equivalently, as 

00 

+ E + + 7nCl] V) = 0. (B.6) 



n=2 



If v is such that a}~_\%v-\ — 0, Eq. (IB.6|) is obviously satisfied by the recurrence relations 
(I117p . If a l _i2/ v -i ^ 0, the above equation is also fulfiled when and or 

and are linearly dependent: these cases give the recurrence relations (11181) and (11191) . 
respectively, as we will see. 

In effect, from (J2"0j) . in a^ 1 ^,_ 1 we have 



ail = ^ W2; o [f + 1 — ^r] 



2() 



B2 

2 



(B.7) 



If there is no cancellation between numerators and denominators, we find a 1 _ 1 = [ol\y^ v -\ 
0) by choosing 



v = v x = (£ 2 /2) - 1 or 1/ = ^ 2 = (Bi/zo) + (5 2 /2). 



(B.8) 



Let us consider only v — V\ corresponding to the solutions Ui given in Eq. f 1 1 1 3 [) (for v = v 2 , 
we get U 2 = 7\Ui). Then, we find that 



V — V\ 



B 2 

— - 1 

2 



0, if B 2 + 1, 2; 
< iw2 [1 + (Bi/zo)] , \IB 2 = 1; 
-iu;z [1 + (£i/*b)] , if B 2 = 2, 



(B.9) 



whereof we see that, in general, oc_]^ v -\ ^ if -B 2 = 1 or B 2 = 2 - for these cases there 
are cancellations between numerators and denominators in (1B.7I) . Thus, we find that Eq. 
( IB. 61) is satisfied by the recurrence relations (11171) for b x n if -B 2 7^ 1,2 (6* = 5* for one-sided 
solutions). Next we will find the relations fl 1 1 8 [) if .B 2 = 1, and (11191) if B 2 = 2. In the first 
place, for the Coulomb functions flA.lOj) we have 



i>n+v{v,y) 



v=(B 2 /2)-l 



v=(B 2 /2)-l 



[2iy] n +^<f> [n + f + i V , 2n + 5 2 ; -2iy] , 
nnlTsljZ^ 1-2WT + ^ [n + f ± 117, 2n + B 2 ; T 2zy] 



where y = uz. Using the relations (m = 1, 2, 3, • • • ) 



$(a,c; j/) 



T(a + m) 



lim 

c _>!_ m r(c) r(a)r(m + i, 

\I/(a, 1 — m; y) = y m \l/(a + m, 1 + m, y) 



y m $(a + m, 1 + m; y) and 



(B.10) 



51 



we find that if B 2 = 1 {y 



-1/2) 



lim • 

m— >— 1 



4 



01, 1p 



± 



4 



± 



u+1- 



Therefore, if B 2 = 1, the first and the third terms in Eq. (IB.6j) are linearly dependent leading 
to the recurrence relations (11181) . In the second place, if B 2 = 2 [y = 0) we find 



lim 

m— >— 1 



-irj(j)o, V-i = 



= -irjfy,, [B 2 = 2, u = 0] 



that is, if B 2 = 2 the first and the second terms in Eq. (1B.6|) are linearly dependent giving 
the recurrence relations (I119p . The previous limits suggest that = if 77 = and, 
consequently, relations f II 1 9 j) reduce to (I117p . In fact, by putting v = and B 2 = 2 in Eqs. 
(p0| we find <2r_i = 0. 



C Other One-Sided Solutions for the CHE 



Now we complete the list of the one-sided solutions for the CHE given in Sec. 3. First we 
write the solutions for n ^ and, by supposing that the series summation begins at n = 
we give the conditions which assures the independence of the Coulomb wave functions as 
well as the recurrence relations for the series coefficients. For series beginning at n > 0, see 
f)128al) and (1128b[) . After this, the solutions for the case rj = are expressed by series of 
Bessel functions. 



C.l Solutions for f] ^ 

Second and sixth sets: 



Conditions: B 2 + ^ ^ -2, -3, -4 • • • ; 



Eqs. (HHDif f + f ^-±,0; (ffH if f + f 



2' 



M) if f + f = 0; 



Uf{z) 



e iulz z 1+ ' 



°E ft 

n=0 



2n +2+B 2 + ^ 



$ |n + zr/ + l + ^f + ^,2n + 2 + fi 2 + ^;-2za;z 



2() 



^1 ^ »n \-2iu}z\' 1 

n=0 

^ [n ± Z77 + 1 + f + ^, In + 2 + 5 2 + t2z'u;z 



(c.r 
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oo 

tulz i+£i bj [2iu(z-z )Y 

U 6 (z) = e lu]Z z *o ^ r ^^;^^l x 



2 n+2+B 2 + ^ 
n=0 L 



$ [n + if] + 1 + + f£, 2n + 2 + 5 2 + -2iu(z - z ) 

oo 

7-°r±/ \ ±iujz 1+^ b n [-2io;(z-*o)]" 

0*(z) = ^ FT^I^Tf] x 

n=0 

* [n ± 27/ + 1 + f + 2n + 2 + 5 2 + t2z'u;(z - z ) 



2<g,z [n+l][n+B a + ^ J 
[2n+2+B 2 + ^Lj [2n+3+B 2 +^Lj ' ^™ Pn ' 

o 6 2^ [n+B 2 + ^L] [n+l+^] [n+tq+gl + ^] H^+j 

Third and seventh sets: 

Conditions: B 2 ^ 4, 5, 6, • • • . 

Eqs. (HZD if 5 2 ^ 2, 3; (HHD if B 2 = 3; dm]) if B 2 = 2. 



oo 
n=0 

$ [n + 27/ + 2 - ff , 2n + 4 - B 2 ; -2iuz] 

oo 

St , „ Bi 



* [n ± ir7 + 2 - ^a, 2n + 4 - B 2 ; T^iouz] 



„ 2ia;z [n+l][n+2+^l] 

rv = - 

(2n+4-B 2 )(2n+5-B 2 ) ' 



Pn = ~{n + l)(n + 2 - fl 2 ) - £3 " 7/^0 - (ir+2-B 2 K2n+ + 4-B 2 J ) 

o 3 2iwz ["+2-B 2 ] [n+l-Ba-fl] [n+l-^+i??] [n+l-^-ir?] 

= (2n+l-B 2 )(2n+2-B 2 ) ' 



(C.2) 



(C.3) 



T[2n+4-B 2 ] 

n=0 

$ [n + ir] + 2 - 2n + 4 - 5 2 ; -2zu;(z - z )] , 

oo (^-'•^0 



n=0 

* [ra ± ir] + 2 - 2n + 4 - B 2 ; ^2iu(z - z )] 



r[n=Fw?+2-4p] 



2^ [n + l][n+2-B 2 -^] » .g 

u ra ( 2n+ 4_ j B 2 )(2n+5-B 2 ) ' ^n> 

o 7 2iwz [n+2-B 2 ] [1+H-f^] [n+l-^+i?;] [n+l-^-M?] 
= (2n+l-B 2 )(2n+2-B 2 ) 
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Fourth and eighth sets: 

Conditions: B 2 + ^ 4 2, 3, 4 ■ • • . 
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C.2 Solutions for 77 = 

The following solutions for the CHE with 77 
one-sided solutions for the spheroidal equation. 
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